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Ich weijJ ehrlich nicht, was die Leute meinen, wenn sie von der Freiheit des men- 
schlichen Willens sprechen. Ich habe zura Beispiel das Gefuhl, dafi ich irgend etwas 
will; aber was das mit Freiheit zu tun hat, kann ich iiberhaupt nicht verstehen. Ich 
spiire, dafi ich meine Pfeife anziinden will und tue das auch; aber wie kann ich das 
mit der Idee der Freiheit verbinden? Was liegt hinter dem Willensakt, dafi ich meine 
Pfeife anziinden 'will? Ein anderer Willensakt? Schopenhauer hat einmal gesagt: 
"Der Mensch kann tun was er will; er kann aber nicht wollen was er will. " 

[Honestly I cannot understand what people mean when they talk about the freedom of 
the human will. I have a feeling, for instance, that I will something or other; but what 
relation this has with freedom I cannot understand at all. I feel that I will to light my 

pipe and I do it; but how can I connect that up with the idea of freedom? What is 
behind the act of willing to light the pipe? Another act of willing? Schopenhauer once 
said: "Man can do what he wills but he cannot will what he wills. "] 



Albert Einstein 
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Abstract 



The descriptions of the quantum realm and the macroscopic classical world differ significantly not 
only in their mathematical formulations but also in their foundational concepts and philosophical 
consequences. The assumptions of a genuine classical world — local realism and macroscopic real- 
ism (macrorealism) — are at variance with quantum mechanical predictions as characterized by the 
violation of the Bell and the Leggett-Garg inequality, respectively. When and how physical systems 
stop to behave quantumly and begin to behave classically is still heavily debated in the physics 
community and subject to theoretical and experimental research. 

The first chapter of this dissertation puts forward a novel approach to the quantum-to-classical 
transition fully within quantum theory and conceptually different from already existing models. It 
neither needs to refer to the uncontrollable environment of a system (decoherence) nor to change 
the quantum laws itself (collapse models), but puts the stress on the limits of observability of 
quantum phenomena due to imprecisions of our measurement apparatuses. Naively, one would say 
that the predictions of quantum mechanics reduce to those of classical physics merely by going to 
large quantum numbers. Using a quantum spin as a model object, we first demonstrate that for 
unrestricted measurement accuracy the system's time evolution cannot be described classically and 
is conflict with macrorealism throTigh violation of the Leggett-Garg inequality, no matter how large 
the spin is. How then does the classical world arise? Under realistic conditions in every-day life, 
we are only able to perform coarse-grained measurements and do not resolve individual quantum 
levels of the macroscopic system. We show that for some "classical" Hamiltonians it is this mere 
restriction to fuzzy measurements which is sufficient to see the natural emergence of macrorealism 
and the classical Newtonian laws out of the full quantum formalism. This resolves the apparent 
impossibility of how classical realism and deterministic laws can emerge out of fundamentally 
random quantum events. 

In the second chapter, we find that the restriction of coarse-grained measurements usually allows 
to describe the time evolution of any quantum spin state by a time evolution of a statistical mixture. 
However, we demonstrate that there exist "non-classical" Hamiltonians for which the time evolution 
of this mixture cannot be understood classically, leading to a violation of macrorealism. We derive 
the necessary condition for these non-classical time evolutions and illustrate it with the example of 
an oscillating Schrodinger cat-like state. Constant interaction of the system with an environment 
establishes macrorealism but cannot account for a continuous spatiotemporal description of the 
system's non-classical time evolution in terms of classical laws of motion. We argue that non- 
classical Hamiltonians are unlikely to appear in nature because they require interactions between 
a large number of particles or are of high computational complexity. 

The third chapter investigates entanglement between collective operators in two specific physical 
systems, namely in a linear chain of harmonic oscillators and in ensembles of spin-^ particles. We 
show that under certain conditions entanglement between macroscopic observables can persist for 
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large system sizes. However, since this analysis uses sharp measurements, it is not in disagreement 
with our quantum-to-classical approach. 

The last chapter addresses the question of the origin of quantum randomness and proposes a 
link with mathematical undecidability. We demonstrate that the states of elementary quantum 
systems are capable of encoding a set of mathematical axioms. Qiiantum measurements reveal 
whether a given proposition is dccidable or undccidablc within this set. We theoretically find 
and experimentally confirm that whenever a mathematical proposition is undecidable within the 
axiomatic set encoded in the state, the measurement associated to the proposition has random 
outcomes. This supports the view that quantum randomness is irreducible and a manifestation of 
mathematical undecidability. 
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Zusammenfassung 



Die Beschreibungen der quantenmechanischen und der klassischen Welt unterscheiden sich nicht 
nur signifikant in ihren mathematischen Formulierungen sondern auch in ihren grundsatzlichen 
Konzepten und philosophischen Implikationen. Die Annalimen einer klassischen Welt - lokaler und 
makroskopischer Realismus (Makrorealismus) - widersprechen den Vorhersagen der Quantenphysik, 
was durch die Verletzung der Bell- und der Leggett-Garg-Ungleichung charakterisiert wird. Die 
Frage, wann und wie physikalische Systeme aufhoren sich quantenmechanisch und anfangen sich 
klassisch zu verhalten, wird in der wissenschaftlichen Gemeinschaft noch immer heftig diskutiert 
und ist Gegenstand experimenteller und theoretischer Forschung. 

Das erste Kapitel der vorliegenden Dissertation entwickelt einen neuen Zugang zum Ubergang 
der Quanten- zur klassischen Physik, und zwar voUkommen innerhalb der Quantentheorie und kon- 
zeptioncU vcrschieden von bcrcits bcstchcndcn ModcUcn. Dicscr Zugang muss sich weder auf die 
unkontroUicrbarc Unigcbung von Systcmcn bczichcn (Dckoharcnz) noch die Gcsctzc der Quan- 
tenmechanik selbst abandern (KoUaps-Modelle). Er fokussiert sich vielmehr auf die Limitierung 
der Beobachtbarkeit von Quantenphanomen aufgrund der Ungenauigkeit unserer Messapparate. 
Naiverweise wiirde man annehmen, dass sich die Vorhersagen der Quantenmechanik auf jene der 
klassischen Physik allein dadurch reduzicrcn, indem man zu groBen Quantenzahlcn gcht. Wir vcr- 
wcnden einen Quantenspin als ModcUobjckt und zcigcn zunachst, dass bci uncingcschranktcr Mess- 
genauigkeit die Zcitcvolution dcs Systems nicht klassisch verstanden werden kann und aufgrund der 
Verletzung der Leggctt-Garg-Ungleichung im Widerspruch zu Makrorealismus steht, selbst wenn 
der Spin beliebig grofi ist. Wie entsteht dann die klassische Welt? Unter realistischen alltaglichen 
Bedingungen sind wir nur in der Lage, grobkornige Messungen durchzufiihren, die die einzelnen 
Quantenniveaus dcs makroskopischen Systems nicht aufloscn konnen. Wir zeigen, dass fiir be- 
stimmtc "klassische" Hamilton-Opcratoren dicsc bloBc Einschrankung zu unscharfen Messungen 
ausreicht, um die natiirliche Emergenz von Makrorealismus und klassischer Newtonscher Gesetze 
aus dem voUen quantenmechanischen Formalismus zu sehen. Dies lost die scheinbare Unmoglichkeit 
auf, wie klassischer Realismus und deterministische Gesetze aus fundamental zufalligen Quantener- 
eignissen entstehen konnen. 

Im zweiten Kapitel zeigen wir, dass die Einschrankung grobkorniger Messungen es iiblicher- 

weise crlaTibt, die zcitlichc Evolution jcdcs bclicbigen quantenmechanischen Spinzustands durch 
die zcitlichc Evolution cincr statistischcn Mischung zu beschreiben. Ungcachtct dcssen demon- 
strieren wir, dass es "nicht-klassische" Hamilton-Operatoren gibt, fiir die die Zeitentwicklung die- 
ser Mischung nicht klassisch verstanden werden kann. Wir leiten die allgemeine Bedingung fiir 
solche nicht-klassischen Zeitentwicklungen her und veranschaulichen sie anhand des Beispiels ei- 
ner oszillierenden Schrodinger-Katze. Andauernde Interaktion des Systems mit einer Umgebung 
etablicrt Makrorealismus, kann aber keine kontimiicrliche raumzcitlichc Bcschreibung der nicht- 
klassischen Zeitevolution des Systems durch klassische Bewegungsgleichungen liefern. Wir argumen- 
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tieren, dass nicht-klassische Hamilton-Operatoren in der Natur wahrscheinlich nicht vorkommen, 
weil sie Vielteilchen-Wechselwirkungen benotigen oder von hoher Komplexitat sind. 

Das dritte Kapitel untersucht Verschrankung zwischen koUektiven Operatoren in zwei spezi- 
fischen physikalischcn Systemen, namlich in einer linearen Kette von harmonischen Oszillatoren 

Tind in Ensembles von Spin-|-Tcilehcn. Wir zcigcn. dass Verschrankung zwischen makroskopischen 
Observablen fiir grol3e Systeme bestehen bleiben kann. Zumal diese Analyse scharfe Messungen 
verwendet, steht sie nicht im Widerspruch zu unserem Zugang zum Ubergang von der Quanten- 
zur klassischen Physik. 

Das letzte Kapitel beschaftigt sich mit der Prage nach dem Ursprung von quant enmechanischem 
Zufall und schlagt eine Verbindung mit mathematischer Unentscheidbarkeit vor. Wir demonstrie- 
ren, dass Zustande elementarer QTiantensysteme einen Satz von mathematischen Axiomen kodieren 
konnen. Quantenmechanische Messungen bringen zum Vorschein, ob eine gegebene Proposition in- 
nerhalb dieses Satzes entscheidbar oder unentscheidbar ist. Wir finden theoretisch und bestatigen 
experimentell, dass Messungen, die mit innerhalb des vom Quantenzustand kodierten Axiomen- 
satzes unentscheidbaren mathematischen Propositionen assoziiert sind, zu zufaUigen Resultaten 
fiihren. Dies unterstiitzt die Sichtweise, dass quantenmechanischer Zufall irreduzibel und eine Ma- 
nifestation von mathematischer Unentscheidbarkeit ist. 
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Chapter 

Introduction 



Since the birth of quantum theory in the 1920s, quantum entanglement and quantum superposi- 
tion have been used to highhght a number of counter-intuitive phenomena. They he at the heart 
of the Einstein-Podolsky-Rosen [30 and the Schrodinger cat paradox (Figure [l]). The corre- 
sponding conflicts between quantum mechanics on one side and a classical world — local realism and 
macroscopic realism (macrorealism) — on the other, are quantitatively expressed by the violation of 
Bell's [9] and the Leggett-Garg inequality [60] , respectively. The quantum violation of local realism 
shows that the view is untenable that space-like separated events do not influence each other and 
that objects have their properties prior to and independent of measurement. The quantum violation 
of macrorealism means that it is wrong to believe that a macroscopic object has deflnite properties 
at any time and that it can be measured without effecting them or their subsequent dynamics. 

The importance of this incongruousness today exceeds the realm of the foundations of quantum 
physics and has become an important conceptual tool for developing new quantum information tech- 
nology. Entanglement and superposition allow to perform certain computation and communication 
tasks such as quantum cryptography [10l|38], teleportation [HJITlj or quantum computation [28|l67j. 
which are not possible classically. Experiments in the near future will be realized with increasingly 
complex objects, either by entangling more and more systems with each other, or by entangling sys- 
tems with a very large number of degrees of freedom. Eventually, all these developments will push 
the realm of quantum physics well into the macroscopic world. Moreover, implications on society in 
a cultural sense may manifest themselves, for the characteristics and peculiarities of the quantum 
world — in particular quantum entanglement and quantum superposition — could eventually become 
part of the every-day experience. 

However, the macroscopic classical world that we perceive around us does not show any char- 
acteristics of the quantum realm. The question "Why do classical systems stop to show quantum 
features?" is still answered in radically different ways within the physics community. Since clas- 
sical apparatuses are needed for performing measurements on quantum systems, this question is 
also related to the so called "measurement problem" and the various interpretations of quantum 
mechanics, ranging from the Copenhagen over the Bohmian to the many- worlds interpretation |46j . 

On the one hand, there exist a number of so-called collapse models fSSlES] which try to explain 
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Introduction 




Figure 1: In 1935, Schrodinger put forward the "burlesque" Gedankenexperiment of a "hell machine" to 
illustrate that, according to quantum physics, it is possible to prepare a cat in a superposition of 'dead' 
and 'alive' [83]. [Picture taken from Sciences et Avenir Hors-Serie 148, 54 (2006).] 

the discrepancy between the quantum and the classical world by introducing a fundamental break- 
down of quantum superpositions at some quantum-classical border. On the other, the decoherence 
program |103l I104j demonstrates that the states of complex systems interacting with an environ- 
ment, which cannot be accessed and controlled in detail, rapidly evolve into statistical mixtures 
and lose their quantum character. 

While neither of these approaches can give a definite or already experimentally settled answer, 
the understanding of the quantum-to-classical transition is not only of prior importance for the 
future development towards macroscopic superpositions and entanglement but also necessary for a 
consistent description of the physical world. Collapse models make assumptions about inherently 
non-quantum mechanical background fields or gravitational mechanisms which are still to be tested 
experimentally. The decoherence program is inherently quantum mechanical and can give good 
explanations for many observations though it has to rely on the assumption of a preferred pointer 
basis. However, the effects of decoherence can in principle be always reduced and the experimental 
progress of the last years has already demonstrated quantum interference of (Schrodinger-cat like) 
macroscopic superpositions, e.g., interference fringes with large molecules of ~ 10'^ atomic mass 
units ID, entanglement between clouds of ~ 10^^ atoms |17], or superpositions of macroscopically 
distinct flux states in superconducting rings corresponding to ~ 10^^ amperes of current flowing 
clock- or anticlockwise [32j . These experiments circumvent the problem of decoherence but did not 
yet come into the region where they could exclude collapse models. 

Until today there exists no definite answer to the problem of the quantum-to-classical transition. 
Hence, certainly one of the most fundamental and interesting questions in modern physics still 
remains unanswered: 

How does the classical physical world emerge out of the quantum realm? 

Chapter 1 of this dissertation addresses this question from a novel perspective and develops 
an approach to the quantum-to-classical transition fully within quantum theory and conceptually 
different from already existing models. It neither needs to refer to the environment of a system 
(decoherence) nor to change the quantum laws itself (collapse models) but puts the stress on 
the limits of observability of quantum phenomena due to our measurement apparatuses. Using a 
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Figure 2: Under the magnifying glass of sharp measurements Albert Einstein sees a strange and colorful 
quantum picture of the face next to him. Its abstractness is symbolized by Pablo Picasso's "Head of a 
Reading Woman" . Under an every-day coarse-grained view the classical appearance of Charlie Chaplin 
emerges. 

quantum spin as a model object, we first demonstrate that for unrestricted measurement accuracy 
the system's time evolution cannot be described classically and is in conflict with macrorealism 
through violation of the Leggett-Garg inequality. This conflict remains even if the spin is arbitrarily 
large and macroscopic. 

Under realistic conditions in every-day life, however, we are only able to perform coarse-grained 
measurements and do not resolve individual quantum levels of the macroscopic system. As we show, 
it is this mere restriction to fuzzy measurements which is sufficient to see the natural emergence of 
macrorealism and the classical Newtonian laws out of the full quantum formalism: the system's time 
evolution governed by the Schrodinger equation and the state projection induced by measurements. 
This resolves the apparent impossibility of how classical realism and deterministic laws can emerge 
out of fundamentally random quantum events. Figure |2] presents an illustration of this approach. 

Chapter 2 flrst shows that a violation of the Leggett-Garg inequality itself is possible for arbi- 
trary Hamiltonians given the ability to perform sharp quantum measurements. Apparatus decoher- 
ence or the restriction of coarse-grained measurements usually allow to describe the time evolution 
of any quantum spin state by a time evolution of a statistical mixture. However, we demonstrate 
that there are "non- classical" Hamiltonians for which the time evolution of this mixture cannot be 
understood classically, leading to a violation of macrorealism. We derive the necessary condition for 
these non-classical time evolutions and illustrate it with the example of an oscillating Schrodinger 
cat-like state. System decoherence, i.e. the continuous monitoring of the system by an environ- 
ment, leads to macrorealism but a dynamical description of non-classical time evolutions in terms 
of classical laws of motion remains impossible. 

In the last part we argue that non-classical Hamiltonians either require interactions between a 
large number of particles or are of high computational complexity. This might be understood as the 
reason why they are unlikely to appear in nature. 

Chapter 3 investigates entanglement between collective operators in two speciflc physical sys- 
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terns, namely in a linear chain of harmonic oscillators and in ensembles of spin-i particles. We 
demonstrate that under certain conditions entanglement between macroscopic observables can in- 
deed persist for large system sizes. However, since this analysis uses sharp measurements, it is not 
in disagreement with our quantum-to-classical approach. 

Chapter 4 addresses the question of the origin of quantum randomness. In our view, classical 
physics emerges out of the quantum world but the randomness in the classical mixture is still 
irreducible and of quantum nature. We propose to link quantum randomness with mathematical 
undecidability in the sense of Chaitin's version of Godel's theorem. It states that given a set of 

axioms that contains a certain amount of information, it is impossible to deduce the truth value of 
a proposition which, together with the axioms, contains more information than the set of axioms 
itself. 

First, we demonstrate that the states of elementary quantum systems are capable of encoding 
mathematical axioms. Quantum mechanics imposes an upper limit on how much information can 
be carried by a quantum state, thus limiting the information content of the set of axioms. Then, we 
show that quantum measurements are capable of revealing whether a given proposition is decidable 
or undccidablc within this set. This allows for an experimental test of mathematical undecidability 
by realizing in the laboratory the actual quantum states and operations required. We theoretically 
find and experimentally confirm that whenever a mathematical proposition is undecidable within 
the system of axioms encoded in the state, the measurement associated to the proposition gives 
random outcomes. Our results support the view that quantum randomness is irreducible and a 
manifestation of mathematical undecidability. 



Chapter 1 



quantum 



Classical 



world emerging from 
physics 



Summary: 



Inspired by the thoughts of Peres on the classical limit |74| — we present a novel theoretical ap- 
proach to macroscopic realism and classical physics within quantum theory. While our approach 
is not at variance with the decoherence program |103} 1104] , it differs conceptually from it. It is 
not dynamical and puts the stress on the limits of observability of quantum effects of macroscopic 
objects, i.e. on the required precision of our measurement apparatuses such that quantum phenom- 
ena can still be observed. The term "macroscopic" is used here to denote a system with a high 
dimensionality rather than a low-dimensional system with a large parameter such as mass or size. 
Furthermore, there is no need to change the quantum laws itself like in collapse models [551 173j . 

Using a quantum spin as a model system, we first show that, if consecutive eigenvalues of 
a spin component can be experimentally resolved in sharp quantum measurements, the Leggett- 
Garg inequality is violated for arbitrary spin lengths and the violation persists even in the limit 
of infinitely large spins. This contradicts the naive assumption that the predictions of quantum 
mechanics reduce to those of classical physics merely due to the fact that a system becomes "large" . 
For local realism this persistence of quantum features was demonstrated by Garg and Mermin [33], 
and the violation even increases with the systems' dimensionality [HJUS]. 

In every-day life, however, one not only encounters very high-dimensional systems but is ex- 
perimentally restricted to coarse-grained measurements. They only distinguish between eigenvalues 
which are separated by much more than the intrinsic quantum uncertainty. We show for arbi- 
trary spin states that, given a certain time evolution, the macroscopically distinct outcomes obey 
the classical Newtonian laws which emerge out of the Schrodinger equation and the projection 
postulate. 

This suggests that classical physics can be seen as implied by quantum mechanics under the 
restriction of fuzzy measurements and resolves the apparent impossibility of how classical realism 
and deterministic laws can emerge out of fundamentally random quantum events. 
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1 Classical world emerging from quantum physics 



This chapter mainly bases on and also uses parts of Reference [56] : 

• J. Kofler and C. Brukner 

Classical world arising out of quantum physics under the restriction of coarse-grained mea- 
surements 

Phys. Rev. Lett. 99, 180403 (2007). 



1.1 The Leggett-Garg inequality 
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1.1 The Leggett-Garg inequality 

In this section we introduce the concept of macroscopic realism (macroreahsm) and show how to 
derive a Leggett-Garg inequahty, which can be used as a tool to indicate whether or not a system's 
time evolution can be understood in classical terms. In agreement with this, we then briefly 
demonstrate explicitly that the inequality can always be violated for genuine quantum systems but 
always satisfied for classical objects. 

Macrorealism is defined by the conjunction of the following three postulates [55]: 

(1) Macrorealism, per se. A macroscopic object which has available to it two or more macroscop- 
ically distinct states is at any given time in a definite one of those states. 

(2) Non-invasive measurahility. It is possible in principle to determine which of these states the 
system is in without any effect on the state itself or on the subsequent system dynamics. 

(3) Induction. The properties of ensembles are determined exclusively by initial conditions (and 
in particular not by final conditions). 

The last two postulates can be phrased into the single assumption that the object's state is 
independent of past and future measurements [SI]. Classical (Newtonian) physics belongs to the 
class of macrorealistic theories. 

Now consider a macroscopic physical system and a dichotomic quantity A, which whenever 
measured is found to take one of the values ±1 only. Further consider a series of runs starting from 
identical initial conditions at time t ~ Q such that on the first set of runs A is measured only at 
times ti and t2, only at ^2 and on the second, at <3 and <4 on the third, and at ti and ^4 on the 
fourth (0 < < ^2 < ^3 < ^a)- Let Ai denote the value of A at time ti (see Figure [TTT]) . Consider 
the algebraic combination of the Clauser-Horne-Shimony-Holt (CHSH) type [23] : 

(A2 - A4) + A3 (A2 + A4) = ±2 . (1.1) 

It can only have the values +2 or —2 for one of the two brackets has to vanish and the other is 
+2 or —2 and then multiplied with +1 or —1. Macrorealism per se is refiected by the objective 
existence of unambiguous values of Ai at all times, and non-invasive measurahility together with 
induction is refiected by the fact that the AiS are the same, independent in which combination 
they appear. E.g., A4 is independent of previous measurements, i.e. whether it appears in a run 

? = ?2 h U 

I 1 1 1 1 ► / 

^1 ^2 ^3 ^4 

Figure 1.1: In a macrorealistic theory macrorealism per se implies that a time-dependent quantity A{t) 
of a macroscopic object has a well defined unambiguous value Ai = A{ti) at every time ti. Non-invasive 
measurahility (together with induction) is reflected by the fact that the Ai^s do not depend on whether the 
system was or was not measured at earlier times. 
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together with Ai or A3. Repeating the experimental runs many times, we introduce the temporal 
correlation functions 

aj = {A,Aj)- (1.2) 



By averaging (1.1) it follows that any macrorealistic theory has to satisfy the Leggett-Garg inequal- 
ity [BD] 

K = Ci2 + C23 + C3i-Cu<2. I (1.3) 



Its violation implies that the object's time evolution cannot be understood classically. 

Let us briefly analyze the quantum evolution of a microscopic quantum object, say the precession 
of a spin-i particle with the Hamiltonian H = ^ uj a^, where ut is the angular precession frequency 
and ax is the Pauli x-matrixj^ If we measure the spin along the z-direction, then we obtain the 
temporal correlations Cij = {uz{ti) az(tj)) = cos[LLi{tj—ti)]. Choosing the four possible measurement 
times as equidistant, with time distance At = t2—ti = ta — ^2 — t^ — ts, the Leggett-Garg inequality 
becomes 

K = 3cos(wAt) - cos(3wAi) < 2 . (1.4) 
This is maximally violated for the time distance At — ^ for which K = 2 \/2 (see red line in 



Figure 1.2 1. The violation is not surprising as a spin-1 particle is a genuine quantum system and 
cannot have the objective properties tentatively attributed to macroscopic objects prior to and 
independent of measurements. In contrast, we consider an arbitrarily sized uniformly rotating 
classical spin vector, again precessing around x and pointing along z at time ti. As dichotomic 
observable quantity we use A{ti) = sgn(coS(jjii) such that A = +1 (—1) if the spin is pointing 



upwards (downwards) along z. As expected, the inequality (1.3) is always satisfied (see blue line in 



Figure 1.2) 



1.2 Violation of the Leggett-Garg inequality for arbitrarily 
large spins 



In this section, we demonstrate that the Leggett-Garg inequality (1.3) is violated for arbitrarily 
large (macroscopic) spin lengths j as long as accurate measurements can be performed. In any run, 
the first of the two measurements acts as a preparation of the state for the subsequent measurement. 
Therefore, the initial state of the spin is not decisive and it is sufficient for us to consider as initial 
state the maximally mixed one: 

p(0)=-^y |m)(TO| = -^. (1.5) 

Here, 1 is the identity operator and |m) are the Jz (spin z-component) eigenstates with the possible 
eigenvalues m — — j, — j + 1, + j. We consider the Hamiltonian 

= ^ + (1-6) 



^Throughout this and the subsequent chapter we use units in which the reduced Planck constant is h= 1. 



1.2 Violation of the Leggett-Garg inequality for arbitrarily large spins 
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Figure 1.2: Violation of the Leggett-Garg inequality |l.4[ ) for a rotating spin-| particle with precession 
frequency lu. The left-hand side of the inequality, K, is shown by a red line, while the classical limit, K — 2, 
is indicated by a dashed line. If the distance At between the four possible equidistant measurement times 
is chosen as At = ^ , then the inequality is maximally violated with K = 2 \/2. The blue line shows the 
left-hand side of the Leggett-Garg inequality for a classical rotating spin vector. Its time evolution can be 
understood classically and does not violate the inequality. 



where J is the rotor's total spin vector operator, its a;-component, / the moment of inertia and 
to the angular precession frequency. The constant of motion ^ can be ignored since commutes 
with the individual spin components and does not contribute to their time evolution. The solution 
of the Schrodinger equation produces a rotation about the x-axis, represented by the time evolution 
operator 



Ut 



(1.7) 



We assume that individual eigenstates |m) can be experimentally resolved and use the parity 
measurement 



^-E (-1) 



j-m 



m) (m| 



(1.8) 



with the possible dichotomic outcomes ±1 (identifying ± — ±1). The temporal correlation function 
between results of the parity measurement A at different (arbitrary) times ti and tj (tj > ti) is 



Cij = Pi+ q]+\i+ + Pi- qj-\i- - Pi+ <l]-\i+ ~ Pi- Qj+li- 1 



(1.9) 



where (Pi-) is the probability for measuring -I- (— ) at and qji\ik is the probability for measuring 
I at tj given that k was measured at ti (k, I = -|-, — ). Furthermore, 



K+ = l-K_ = i((itJ-f 1), 
9j + |i:± = 1 - 1]-\i± = \ ((^tj>± + 1) ■ 



(1.10) 

(1.11) 



Here {At- ) is the expectation value of A at ti and {At- ) ± is the expectation value of A at tj given 
that ± was the outcome at ti. The totally mixed state is not changed until the first measurement: 
p{t.i) = ilt^ /5(0) Ul /5(0) and we find 



{At^)^Tr[mA] 



2j 



0. 



(1.12) 
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The approximate sign is accurate for half integer j as well as in the macroscopic limit j 3> 1, which 
is assumed from now on. Hence, we have pi+ = pi- = \, which is self-evident for a totally mixed 
state. Depending on the measurement result ± at t^, the state is reduced to 

P.(t.) = ^^^^ = ^, (1.13) 
Tr[P±p(t,)P±] 2j + l' ^ ' 

with P± = ^ (1 ± A) the projection operator onto positive (negative) parity states. Denoting 
9 = UJ {tj — ti), the remaining expectation value {At.)± — Tr[iJt._t. p±{ti) tj}^_^. A] becomes 

(At.)± = ±— Tr[e-'^-^"=^ e'^'^J-J^) ^^^U-J.)^ 

^ ^ Trje^''^-^"-] ^ ^ sin[(2j + 1) u; (t, - t,)] 
2j + l {2j + l)sm[Lo{t,-t{)] ' 

Here we used the geometrical meaning of the rotations in the first line. From {At.)+ = —{At-)- 
it follows q+\+ + q+\- — 1. Using this and pi+ — \ from above, the temporal correlation function 
becomes 

(1.15) 

Having four possible equidistant measurements with time distance Ai, and using the abbreviation 

a;=(2j + l)a;At (1.16) 



the Leggett-Garg inequality ( |1.3[ ) now reads 

^ nr sin .'^t' 

(1.17) 



3sin.T sinSx 

K w < 2. 

X "ix 



We approximated the sine function in the denominator, assuming 2j+i ^ ^- Inequality (1.171 is 
violated for all positive x < 1.656 and maximally violated for x sa 1.054 where K w 2.481 (compare 
with Reference |74j for the violation of local realism) as can be seen in Figure 1.3 For every spin 
size J, and given a precession frequency tj, it is always possible to choose the time distance At such 
that K >2. 

We can conclude that a violation of the Leggett-Garg inequality is possible for arbitrarily 
high- dimensional systems and even for initially totally mixed states. 

Note, however, that the temporal precision of our measurement apparatuses, which is required 
for seeing the violation, increases with j, as u At has to scale with in order to keep a; « 1. 
Moreover, due to the nature of the parity measurement, consecutive values of m have to be resolved. 



1.3 The quantum-to-classical transition for a spin-coherent 
state 

In this section we will show that coarse-grained measurements not only lead to the validity of 
macrorealism but even to the emergence of classical physics for a certain class of quantum states. 
The generalization to arbitrary states will be done in the next section. 



1.3 The quantum-to-classical transition for a spin-coherent state 
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Figure 1.3: Violation of the Leggett-Garg inequality (1.171 for arbitrarily large spin size j. The left hand 



side of the inequality, K, is indicated by a red line and the classical limit is if = 2 (dashed line). The 



initial quantum state is totally mixed, eq. (1.5 I, and the Hamiltonian (1.61 produces a precession around 



the a;-axis with frequency lo. The dichotomic observable is the parity in a measurement (1.8 1 of the spin's 



z-component. For all spin sizes j and precession frequencies ui, one can always find a time distance At 
between the four possible equidistant measurement times such that the quantity x = (2j + 1) a; At is around 
the value 1 and the classical limit is violated. The maximal violation K ~ 2.481 is achieved for x ~ 1.054. 



Let us start with a preliminary remark about distinguishability of states in quantum theory. 
Any two different eigenvalues mi and TO2 in a measurement of a spin's z-component correspond to 
orthogonal states without any concept of closeness or distance. In Hilbert space the vectors |m) 
and |m+l) are as orthogonal as |m) and |m+10^"): 

(m+l|m)=0, (1.18) 
(m+10i°|m) = 0. (1.19) 

The terms "close" or "distant" only make sense in a classical context, where those 
eigenvalues are treated as close which correspond to neighboring outcomes in the real 
configuration space. 

For example, the "eigenvalue labels" m and m + 1 of a spin component observable correspond 
to neighboring outcomes in a Stern-Gerlach experiment. (Such observables are sometimes called 
classical or reasonable |101lH^F7^ .) If our measurement accuracy is limited, it is those neighboring 
eigenvalues which we conflate to coarse-grained observables. It seems thus unavoidable that certain 
features of classicality have to be assumed beforehand to give the Hilbert space some structure 
which it does not have a priori. 

In what follows, we will first consider the special case of a single spin coherent state and then 
generalize the transition to classicality for arbitrary states. Spin-j coherent states \ = |78}l6| 
are the eigenstates with maximal eigenvalue of a spin operator pointing into the direction = (i?, V') i 
where and tp are the polar and azimuthal angle in spherical coordinates, respectively: 



(1.20) 
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Let us consider the initial spin coherent state at time t = pointing into the direction (t?o, ipo). In 
the basis of Jz eigenstates it reads 

1^0, Vo) = LVmY' cos^'+'" t sin^-™ t e-™'^" |m) . (1.21) 



Under time evolution Ut = e "^'•'==, eq. (1.7 1, the probability that a measurement at some later 



time t has the particular outcome m is given by the binomial distribution 

p{m,t)^\{m\^t,^t)\^ (1.22) 

with cosi3t = sin tjt sin i^o sin (^0 + cos cos z9oi tan 1^94 = cosajitan(/3o ^ sinti;ttan^"'^z9ocos~'"'^(poi 
where "dt f^nd ipt are the polar and azimuthal angle of the (rotated) spin coherent state jiJt, ipt) = 



Ut \'0o,ipo) at time t. In the macroscopic limit, j 3> 1, the binomial distribution (1.221 can be very 
well approximated by a Gaussian distribution 

1 (m-M)^ 

p{m,t)K e (1.23) 

V27r (7 



with a = \fYF^ Bvei'dt the width (standard deviation) and \x = jcos'dt the mean value. 

Under the "magnifying glass" of sharp measurements individual eigenvalues m can be distin- 



guished and the Gaussian probability distribution p{m, t) can be resolved, as shown in Figure 1.4 ^a), 
allowing a violation of the Leggett-Garg inequality. Let us now assume that, as in every-day life, 
the resolution of the measurement apparatus. Am, is restricted and that it subdivides the 2j + 1 
possible different outcomes m into a much smaller number of coarse-grained "slots" m. If the 
slot size is much larger than the standard deviation a ^ ^/J, i.e. 



Am>Vj, (1.24) 



the sharply peaked Gaussian cannot be distinguished anymore from the discrete Kronecker delta. 

Am > y/j : p{m, t) drn,p, ■ (1-25) 

Here, m is numbering the slots (from — j + to j — in steps Am) and ft is the number of the 
slot in which the center /i of the Gaussian lies. This is indicated in Figure [L4jb) . 

In the classical limit, j — s- 00, one can distinguish two cases: (1) If the inaccuracy Am scales 
linearly with j, i.e. Am = 0{j), the discreteness remains. (2) If Am scales slower than j, i.e. 
Am = o(j) but still Am 3> V7, then the slots seem to become infinitely narrow. Pictorially, the 
real space length of the eigenvalue axis, representing the 2j + 1 possible outcomes m, is limited in 
any laboratory, e.g., by the size of the observation screen after a Stern-Gerlach magnet, whereas the 
number of slots grows as j /Am. Then, in the limit j 00, the discrete Kronecker delta becomes 
the Dirac delta function. 

Am ^ & j — > cx) : p{m,t) ^ 6(m — fl) , (1-26) 
which is shown in Figure [L4[c). 



1.3 The quantum-to-classical transition for a spin-coherent state 
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Figure 1.4: An initial spin-j coherent state li^OiV'o) precesses into the coherent state \'&t,ifit) at time t 
under a quantum time evolution, (a) The probability p(m, t) for the outcome m in a measurement of the 
spin's 2-component is given by a Gaussian distribution with width a and mean /i, which can be seen under 
the magnifying glass of sharp measurements, (b) The measurement resolution Am is finite and subdivides 
the 2j + 1 possible outcomes into a smaller number of coarse-grained "slots" . If the measurement accuracy 
is much poorer than the width a, i.e., Am ^ y^, the sharply peaked Gaussian cannot be distinguished 
anymore from the discrete Kronecker delta 5m, p where m is numbering the slots and fi is the slot in which 
the center /i of the Gaussian lies, (c) In the limit j ^ oo the slots seem to become infinitely narrow and 
5m, M becomes the Dirac delta function 5{m — p,). 



Now we have to focus on the question in which sense coarse-grained von Neumann measurements 
disturb the spin coherent state. Let 

,,-J"i)(™l (1-27) 

denote the projector onto the slot fh with {m} the set of all m belonging to m. Then Pm I'd, is 
almost {"djifi) (the zero vector 0) for all coherent states lying inside (outside) the slot, respectively: 

£. iq \^/ 1^'^) for inside m, 
I L) lor /i outside m. 

This means that the reduced (projected) state is essentially the state before the measurement or 
projected away. If ip) is centered well inside the slot, the above relation holds with merely expo- 
nentially small deviation. Only in the cases where j-i?, (p) is close to the border between two slots, the 
measurement is invasive and disturbs the state. Presuming that the measurement times and/or slot 
positions chosen by the observer are statistically independent of the (initial) position of the coher- 
ent state, a significant disturbance happens merely in the fraction a/ Am <C 1 of all measurements. 
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This is equivalent to the already assumed condition -y^y <^ Am. Therefore, fuzzy measurements of a 
spin coherent state are almost always non-invasive such as in any macrorealistic theory, in particular 
classical Newtonian physics. Small errors may accumulate over many measurements and eventually 
there might appear deviations from the classical time evolution. This, however, is unavoidable in 
any explanation of classicality gradually emerging out of quantum theoryj^ 

Hence, at the coarse-grained level the physics of the (quantum) spin system can completely be 
described by a "new" formalism, utilizing an initial (classical) spin vector 3 at time t = 0, pointing 
in the (tJq, ¥'o)-direction with length J = |J| = \/j(i + l) ~ j, where j ^ 1, and a (Hamilton) 
function ^ 

H='^ + ujJ,. (1.29) 

At any time the probability that the spin vector's z-component Jcos'i?t ~ jcosi^t is in slot m is 
given by Sm.fii eq. (1.25), as if the time evolution of the spin components Ji {i = x,y,z) is given 



by the Poisson brackets, 

j^ = [J^,H]pB, (1.30) 



and measurements are non-invasive. Only the term lu Jx in eq. (1.291 governs the time evolution 
and the solutions correspond to a rotation around the x-axis. The spin vector at time t points in the 
(i^t, (^()-direction where iDt and ipt are the same as for the spin coherent state and the probability 



of measurement outcomes is given by S{m — fi), eq. (1.261 



This is classical (Newtonian) mechanics of a single spin emerging from quantum physics. 



1.4 The quantum-to-classical transition for an arbitrary spin 
state 

Now we demonstrate that the time evolution of any spin-j/ quantum state becomes classical under 
the restriction of coarse-grained measurements. At all times any (pure or mixed) spin-j density 
matrix can be written in the quasi-diagonal form [3] 



J I p{n) \n){n\(fn (1.31) 



with d^il = sin I? di? d(y9 the infinitesimal solid angle element and P{il) a not necessarily positive 
real function with the normalization // P{^) = 1- 



The probability for an outcome m in a measurement in the state (1.31 1 is given by 



w{m) = Tr[p |m) (m|] = / / P{n) p{m) d^n , (1.32) 



^For the general trade-oflf between measurement accuracy and state disturbance for the more reaUstic smoothed 
positive operator value measure (POVM) and for related approaches to classicality, see References 1191 1771 135| . A 
natural way of implementing coarse-grained measurements as POVM is presented in the next chapter. In contrast 
to von Neumann measurements, they do not allow to distinguish perfectly between two states at two sides of a slot 
border. But under all circumstances it is in general unavoidable that quantum measurements are invasive to some 
extent. Classicality arises in the sense that the effects of these deviations become negligibly small. 



1.4 The quantum-to-classical transition for an arbitrary spin state 
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Figure 1.5: Given a coarse-grained measurement of the spin's z-component (Am one can compute 

the probabihty Wfn that the outcome wiU be a certain slot m via integration of the quantum state's P- 
function over the region flfn- This is the region between two circles of latitude corresponding to the slot 
m of Jz eigenstates. The coarse-graining condition Am ^ ^ ensures that alternatively the integration 
can be carried out over the state's Q-function, representing the probability distribution of an ensemble of 
classical spins. 



where p(m) is written in eq. (1.221. At the coarse-grained level of classical physics only the proba- 



bility for a slot outcome m can be measured, i.e. 

Tr[pPrn] = y] ^,w{m) 



(1.33) 



with P.„, from eq. (1.271. Inserting eq. (1.31), we get = JJ P{n)Ti[P,n\n){n\]d^n. Using 



eq. (1.281, Am ^Vj, and Tr[|ri)(ri|] = 1, this can be well approximated by 



where is the region between two circles of latitude at polar angles i?i(m) and i?2(w) correspond- 



ing to the slot m (Figure 1.51. We will show that can be obtained from a positive probability 
distribution of classical spin vectors. Consider the well known Q-function [1] [2] 



Qin) = J j Pin') cos^^ f d^n' 



(1.35) 



with d^n' = sin -d' d-d' dip' and 6 = 2arccos{i [1 -I- cos?? cos??' -I- sin??sin?9'cos(v9 - v')]}^^^ the 
angle between the directions = (??, (ys) and fi' = {j)' ,Lp'). In the case of large spins the factor 
cos**^ ® ~ exp(— jG'^) in the integrand is sharply peaked around vanishing relative angle O and 
significant contributions arise only from regions where < l/VJ- The normalization factor "^^^ 



in eq. (1.351 is the inverse size of the solid angle element for which the integrand contributes 



significantly and makes Q normalized: jj^Q{n) d^n — 1. 

The distribution Q is positive because it is, up to a normalization factor, the expectation value 
of the state jfi): 

^ ^ {n\p\n) . (1.36) 



Q{n) = 



An 
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For fuzzy measurements with (angular) inaccuracy A@ ~ 1^2 ("i) — i?i(m) ^ lA/j, which is equiv- 
alent to Am ^ y/j, the probability for having an outcome rh can now be expressed only in terms 
of the positive distribution Q: 

J J Q{n)d^fl. (1.37) 



Figure 1.5 shows the integration region flfn over which P and Q have to be integrated. The 



approximate equivalence of eqs. (1.341 and (1.37) is verified by substituting eq. (1.351 into (1.371 



and is not accurately fulfilled for quantum states p directly at a slot borderj^ 

Note that Q is a mere mathematical tool and not experimentally accessible in coarse-grained 
measurements. Operationally, because of Am ^ y/j an averaged version of Q, denoted as R, is 
used by the experimenter to describe the system. Mathematically, this function R is obtained by 
integrating Q over solid angle elements corresponding to the actual measurement inaccuracy. In the 
classical limit, without the "magnifying glass" , the regions given by the experimenter's resolution 
become "points" on the sphere where R is defined. 

Thus, under coarse-grained measurements, a full description of an arbitrary quantum 
spin state is provided by an ensemble of classical spins with a positive probability distri- 
bution. 



In other words, there exists a hidden variable description. The time evolution of the general 



state (1.31 1 is determined by (1.6 1. In the classical limit it can be described by an ensemble of 



classical spins characterized by the initial distribution Q (i?) , where each spin is rotating according 



to the Hamilton function (1.291. From eq. (1.371 one can see that for the non-invasiveness at the 



classical level it is the change of the distribution Q ( R) which is important and not the change of 
the quantum state or equivalently P itself. In fact, upon a fuzzy measurement the state p is 
reduced to one particular state depending on the outcome rn. 



(1.38) 



with the corresponding (normalized) functions Pm, Qm and The reduction to pm happens 



with probability Wm, which is given by eq. (1.341 or (1.37). Whereas the P-function can change 



dramatically upon reduction, Qm is (up to normalization) approximately the same as the original 
Q in the region ilm. Thus, 



Q,-„(r!) cx {n\p^\n) oc {n\PrnpP,n\^) » {n\p\n) « 



Q(ri) for inside il^, 
for Q, outside fi^- 



(1.39) 



Therefore, at the coarse-grained level the distribution Qm (Rm) of the reduced state after the mea- 
surement can always be understood approximately as a subensemble of the (classical) distribution 
Q before the measurement. 



■^This issue is related to the fact that our von Neumann coarse-grained measurements have sharp slot borders and 
will be resolved in the subsequent chapter. 



1.5 An alternative derivation 
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Effectively, the measurement only reveals already existing properties in the mixture and 
does not alter the subsequent rotation of the individual classical spins. 

The disturbance at the slot borders at that level is quantified by how much Qm differs from a 
function which is (up to normalization) Q within and zero outside. One may think of dividing all 
quantum states and their Q-distributions into two extreme classes: The ones which show narrow 
pronounced regions of size comparable to individual coherent states and the ones which change 
smoothly over regions larger or comparable to the slot size. The former can be highly disturbed 
but in an extremely rare fraction of all measurements. The latter is disturbed in general in a single 
measurement but to very small extent, as the weight — in terms of the Q-distribution — on the slot 
borders (oc-^/J) is small compared to the weight well inside the slot (cx Am). (In the intermediate 
cases one has a trade-off between these two scenarios.) The typical fraction of these weights is 
yj/Am ^1. In any case, classicality arises with overwhelming statistical weight. In the next 
chapter we will introduce measurements with smooth borders (POVM) and therefore reduce the 
disturbance dramatically, even for states near a slot border. 

Finally, we want to point out explicitly: The angular resolution which is necessary to see the 
quantumness, i.e. the superposition character, of a given quantum state is of the order of l/v^- In 
other words, it is necessary to be able to distinguish at least of the order of\/J different measurement 
outcomes. For a macroscopic object, j ^ 10^°, it would be necessary to resolve ~ 10"'^° different 
measurement outcomes. If this precision cannot be met, macrorealism emerges out of quantum 
physics for the rotation Hamiltonian. 



1.5 An alternative derivation 

For the sake of completeness we now present an alternative way to derive that classicality emerges 



under coarse-grained measurements. Again, we consider the totally mixed state (1.5 1 and the 



time evolution (1.7 1. We remind that this allows to violate the Leggett-Garg inequality if sharp 
measurements can be performed. Now we are interested in the probability for obtaining the results 
mi at time ti and m2 at t2 in measurements of the spin operator's z-component — in analogy 
to |74j . where a generalized singlet state and correlations in space are considered. This probability 
can be written as 

p{mi,ti;m2,t2) = p(mi, ti) p(m2, t2)mi,ti , (1-40) 

i.e. as the probability that mi is obtained at ti times the probability that m2 is the result at t2 
given mi at ti. Let = |m)(m| denote the projector onto |m). Using that p{ti) = /5(0), we 
have p{mi,ti) = Tr[/5(ti)PmJ = 2J+1' ^effecting the fact that all of the 2j + 1 eigenstates are 
equally probable in a maximally mixed state. If mi was obtained at ti, the state is reduced to 
Pmi{ti) = Pmi = \'mi){mi\. Then, with At = t2 - ti and 9 = iv At, we have p(m2, t2)mi,ti = 
Tr[Ut,-t,Pm, Ul-t An2] - 1(^2 le-'^-^^ |mi)|2. Therefore, 

p{mi,ti;m2,t2) = — ^ [(msl e""''^- |mi)|2. (1.41) 
2j + 1 
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To continue we perform a discrete Fourier transform with "frequencies" ^ and 77: 
Pit V) = y] e'«™i+''"^) p(mi, ii; m2, h) 

' ^mi ,7712 

Tr[e''>'^- e"'^-^- e'^-^- e'^'^-] Tr[e-^'''-'^] sin[(2j + 1) f ] 



2j + l 2j + l (2j + l)sinf 



(1.42) 



where the four rotations were expressed as a single rotation e~"* '^ around a vector k about an 
angle of the size k — \k\ and the trace was evaluated in the basis where /« • J is diagonal. Following 
firmly Reference [71], it is enough to treat the rotation matrices as 2 x 2 as if j = i, since j does 



not affect the geometrical meaning. If one equates the product of the four matrices in (1.421 to 
g-iK, cr/2 (^-^^rith & the vector of Pauli matrices), then one obtains the dependence of k on ^, 77 and 9: 

cos ^ — cos I cos ^ — sin I sin ^ cos . (1-43) 

Due to noise or coarse-graining of measurements, respectively, the high frequency components 
(^,77 ~ 1) are not experimentally observable. Because of ^,77 ^ high frequency components 
correspond to a sharp measurement resolution Am. The low frequency limit (^,77 <^ 1) of the 
expression for k in the first non-trivial order reads 

«^^ + 77^ + 2^77Cose'. (1.44) 

Now consider a classical spin vector 3 at time t — with length J = |J| = \/ j{j + 1 ) ~ j + \, 
where j ^ \. In analogy to the quantum case, the probability for measuring the classical spin 
vector's z-component as 7771 at ti and 7772 at ^2 is 

Pci.j{rni,ti;m2,t2) = Pc1,j(toi, ti) pcLj(»772, i2)mi,ti • (1-45) 

The first probability is given by Pci,j(wi,ii) = S{mi — • J(ii)), where is the unit vector in 
z-direction and 5 is Dirac's delta function. The classical Newtonian time evolution of the spin 



components Ji {i = x,y,z) is given by the Poisson brackets (1.301, where the (classical) Hamilton 
function is given by (1.291. As 3^ = j (j + 1) is a constant of motion, only the term uj J^. governs 



the time evolution. As we already know, the solutions correspond to a rotation around the a;-axis. 
Thus, the spin at ti reads J(<i) = Rx{i-u ti) 3(0), with Rx{l) the 3x3 rotation matrix about the 
X-axis by an angle 7. But for the scalar product ■3{ti) in Pci,j(''77i, ti) it does not matter whether 
we rotate J(0) around xhy ldIi or around x by —ujti. Hence, we can write 

Pci,j(mi, ii) = 5{mi - a ■ J(0)) , (1.46) 

where a = Rx^—Loti) is the unit vector with polar angle uti and azimuthal angle — ^. The 
remaining probability is 

Pcl,J{m2,t2)r,^^,t^ = (5(7772 - /3 • J(0)) , (1.47) 

i.e. it depends only on the initial spin J(0) and the elapsed time t2 but not on anything that 
happened at ti because of the non-invasiveness of classical measurements. In other words, there 
is no reduction at ti and the condition "777,1 was the outcome at ii" is unnecessary. Here, /3 — 



1.5 An alternative derivation 
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Rx{—LiJ ^2) ©2 is the unit vector with the polar angle lo t2 and azimuthal angle — ^. Hence, eq. ( 1.45 1 
becomes 

Pci,j("^i,ii; m2,t2) = S{mi OL 3) 5{m2 - /3 • J) (1-48) 
with J = J(0). The Fourier transform reads 



with k = ^ a + ry /3. This implies 



C + rf + 2£,r^cos9 , 



(1.49) 



(1.50) 



where 9 = lo {t2 — ti) is the angle between a. and (3, in total agreement with the low frequency limit 
(1.441. The classical correlation Pci,j(^,77) must be averaged over all possible initial directions of 
J = J(f2). Mimicking the mixed quantum state, the distribution is isotropic and thus 



1 

47r 



PcLJ(C,??)d'^^ = 



sin[(2j + l) 



(1.51) 



(2j + l)f 

with d^f2 the infinitesimal solid angle element, the vector k in Pci.3{£,iii) ^ the polar integration 
axis and J ^ This is the limiting value of the quantum correlation (1.421. We note that both 



the quantum and the classical correlation for measurements of a rotating spin in time, eqs. (1.421 
and ( |1.51| ), have a similar form as in the case of a generalized singlet state and measurements in 
space ^74) . 

The quantum to classical transition. Let us now evaluate under which circumstances the quan- 
tum correlation (1.42) becomes classical (1.511, i.e. 

Pel- 



P 



(1.52) 



First, we note that the allowed frequencies ^ and j] are independent and thus have to be small in 
the same order, i.e. ^ ^ ij, whereas higher frequencies are cut off. In the following, for the sake of 
a short and intuitive notation, we consider all quantities as positive and ignore any factors of the 
order of 1 whenever we write ^. Comparing the exact expression for k^, obtained from eq. (1.431, 
with the low frequency classical limit fc^, eq. (1.50), the leading order of the error is 

,2 



Since /c ~ ^ due to eq. (1.50), we have 



k+e 



(1.53) 



(1.54) 



Thus, K — > A: is fulfilled if and only if ^ k, or equivalently, ^ 1. Two formal conditions — (i) for 
the arguments in numerators and (ii) for the denominators in eqs. (1.42) and (1.51), respectively — 
have to be met to guarantee p ^ Pd- 



(2j + 1) 



(2j + 1) 



(2j + l)sinf ^(2j + l)| 



(1.55) 
(1.56) 
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1 Classical world emerging from quantum physics 



Using eq. (1.54), the left-hand side of condition (1.551 becomes 



(2i + i)f -(2j + i)|+je^ 



(1.57) 



As j is very large in the macroscopic limit, we have neglected the term ^"^ compared to j which is 



the leading order of the error. Eq. ( 1.57) is a good approximation for {2j + 1) |, i.e. the right-hand 
side of (1.55), if and only if the error j is much smaller than the smallest term in {2j 



which is I 



^. Hence we have to postulate j <C ^, or equivalently. 



1 

71 



1) 



(1.58) 



and the same condition has to hold for the frequency rj. The evaluation of condition (1.56), using 



a Taylor expansion of the sine function, leads to the same condition (1.58). This high-frequency 
cut-off, ensuring p — > pd, implies that different values of m can be experimentally distinguished 
only if their separation Am fulfills 

Am > (1-59) 

which is the minimum quantum uncertainty for spin coherent states and in agreement with the 
previous sections. 

// the angular resolution of the instruments is much poorer than the intrinsic quantum 
uncertainty, they cannot detect the quantum features of the spin system, let alone a 
violation of the Leggett-Garg inequality. The temporal correlations become classical. 



Chapter 2 



General conditions for quantum 
violation of macroscopic realism 

Summary: 

We first show that a violation of the Leggett-Garg inequality itself is possible for arbitrary Hamil- 
tonians given the ability to distinguish consecutive eigenstates in sharp quantum measurements. 
This is understandable because it is generally accepted that "microscopically distinct states" do not 
have objective existence. For testing macrorealism one needs to apply the Leggett-Garg definition 
referring to macro scopically distinct states. In our every-day life, to experience macrorealism it is 
usually sufficient to employ apparatus decoherence (where the system is isolated and only after a 
premeasurement, i.e. coupling of system and apparatus, the environment interacts irreversibly with 
the apparatus) or the restriction of coarse-grained measurements. 

Both mechanisms usually allow to describe the time evolution of any quantum spin state by a 
classical time evolution of a statistical mixture. However, we demonstrate that there exist "non- 
classical" Hamiltonians for which the time evolution of this mixture cannot be understood classi- 
cally. Despite the fact that apparatus decoherence or coarse-graining allow to describe the state 
merely by a classical mixture at every instance of time, a non-classical Hamiltonian builds up su- 
perpositions of macroscopically distinct states and allows to violate macrorealism. We find the 
necessary condition for these non-classical evolutions and illustrate it with the example of an oscil- 
lating Schrodinger cat-like state. System decoherence, i.e. the constant monitoring of the system by 
an environment, leads to macrorealism but a continuous spatiotemporal description of non-classical 
time evolutions in terms of classical laws of motion remains impossible. 

In the last part we argue that non-classical Hamiltonians require interactions between a large 
number of particles or are computationally much more complex than classical Hamiltonians, which 
might be the reason why they are unlikely to appear in nature. 

This chapter mainly bases on and also uses parts of References [531 ISZ] '■ 
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2 General conditions for quantum violation of macroscopic realism 



• J. Kofler and C. Brukner 

A coarse-grained Schrodinger cat 

In: Quantum Communication and Security, ed. M. Zukowski, S. Kilin, and J. Kowalik (lOS 
Press 2007). 

• J. Kofler and C. Brukner 

The conditions for quantum violation of macroscopic realism 
Phys. Rev. Lett, (accepted); arXiv:0706.0668 [quant-ph]. 



2.1 Violation of the Leggett-Garg inequality for arbitrary Hamiltonians 
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2.1 Violation of the Leggett-Garg inequality for arbitrary 
Hamiltonians 



In contrast to the Leggett-Garg inequality in the previous chapter, ineq. (1.3), with four possi- 
ble measurement times ti, we will now only use three. Then, any macrorealistic theory predicts 
a Leggett-Garg inequality of the Wigner type [55], where Cij = (AiAj) denotes the temporal 
correlation of a dichotomic quantity A at times U and tj : 

K=Ci2 + C23 - Ci3 < 1 . (2.1) 

We extend the approach of Peres in Reference [74] and look at the "survival probability" of the 
system's initial state at time t = 0|^ This state be denoted as |V'(0)) = IV'o) (which must not be an 
energy eigenstate) and, without measurements, it evolves to 

\m) - e-'^* IV-o) (2.2) 

according to the Schrodinger equation in units where the reduced Planck constant is h — 1. Our 
dichotomic observable is 

iEE2|V'o)(V'o|-l, (2.3) 
i.e. we ask whether the system is (still) in the state \tpo) (outcome + = +1) or not (outcome 



-1). The temporal correlations Cij, eq. (1.9l, can be written as 



dj = pi+ qj+\i+ + Pi- qj-\i- - Pi+ q]~\i+ - Pi~ <lj+\i- , (2.4) 

where pi^ (Pi-) is the probability for measuring + (— ) at t.^ and qjmk is the probability for measuring 
I at tj given that k was measured at ti (fc, I = +, For simplicity we choose = and equidistant 



possible measurement times with At = t2 — ti = t^ ~ t2 (Figure 2.1 1. Then the correlation C12 is 
given by C12 = 2p{At) — 1, where 

p{t)^\{yjo\m)\' (2.5) 

is the (survival) probability to find \tpo) given the state \ip{t)). Analogously, we find C13 = 2p{2At) — 
1. As P2+ 1 - P2- = p{At), (73+12+ = 1 - 93-12+ = p(At), and ^3-12- = 1 - 93+12-, the only 
remaining unknown quantity in C23 is g3_|_|2_. For its computation one needs the reduced state at t2, 
given that the outcome was — , \ip-{t2))- It is obtained by applying the projector 1 — |V'o)(V'ol to the 
state at t2, |i/'(i2)), and normahzing: |-0_(i2)) — [|V'(^2))~(V'o|'</'(^2)) \'4'a)]/\/^—p{At). This evolves 
a time At to t^, resulting in Itp^it^)) = Mts)) - (V'o|V'(t2)) \^p{t2))]/x/l-p{At), and 93+12- = 



|('0o|V'-(^3))P is the probability for the outcome + in that state. Plugging everything into (2.11, 
one ends up with the Leggett-Garg inequality 



K ^ 'lp{At)y/p{2At) cos7-4p(2At) + l < 1^ 



(2.6) 



where 7 = 2a — /? and a and /3 are the phases in (V'olV'(^2)) = \/p{^t) exp(iQ;) and {'ipoli^i'ts)) — 
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2 General conditions for quantum violation of macroscopic realism 



At At 

Figure 2.1: At time ti = the initial state of the quantum system is jV'o)- Possible later measurement 
times are t2 — At and = 2At, where the system can be asked whether it is still in its initial state or in 
the orthogonal subspace. 

exp(i/3). 

Now, independent of the system's dimension, it is sufficient to consider as initial state a superpo- 
sition of only two energy eigenstates and \u2) with energy eigenvalues Ei and £"2, respectively: 



IV'o) = (|wi) + |m2))/%/2. Ineq. (2.6l becomes 

K = 2 cos{AEAt) - cos(2A£;A<) < 1 , (2.7) 

with AE = E2 — El the energy difference of the two levels. 

Any non-trivial Hamiltonian leads to a violation of this inequality. 

The left-hand side reaches if = 1.5 for At = and At — ^^i; ^^^d in periods thereof 
(Figure [Q]). 



2.2 Macrorealism per se 

If every Hamiltonian leads to a violation of the Leggett-Garg inequality, why then do we not see 
this in every-day life? The possible answers within quantum theory are: 

1. It is due to decoherence, i.e. the quantum system interacts with an uncontrollable environment 
such that it is driven into a statistical mixture |1031 1104] . 

2. The fact that the resolution of our measurement apparatuses is not sharp, makes it impossible 
to project onto individual states and to see the above demonstrated violation that is always 
present for microstates 

For testing macrorealism — i.e. testing the Leggett-Garg inequality under the restriction of 
coarse-grained measurements — we consider again a spin-j system (with j 3> 1) as a model ex- 
ample. To keep this chapter self-contained, we briefly repeat that any spin-j state can be written 
in the quasi-diagonal form 

JJp{n)\n){n\d^n (2.8) 

^In Reference i74i . exercise 12.23 correctly claims that any non-trivial Hamiltonian is incompatible with a Leggett- 
Garg inequality. But the given reason is inconclusive (and no proof is presented). The survival probability, to which 
is referred to, is written only in form of an inequality, ineq. (12.136) in |74l . This constraint is not strong enough 
and does not exclude time evolutions which would not violate any Leggett-Garg inequality. 



2.2 Macrorealism per se 
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K 




Figure 2.2: Violation of the Wigner-type Leggett-Garg inequality (2.7 1 for arbitrary Hamiltonian evolu- 
tions. Every non-trivial time evolution is in conflict with a classical description as long as one can project 
onto individual quantum states of the system. 



with d^ri the soHd angle element and P a normahzed and not necessarily positive real function [3]. 
The spin coherent states |f2) = with ■& and f the polar and azimuthal angle, are the eigen- 

states with maximal eigenvalue of a spin operator pointing into the direction f2 = |78L 15]: 

Jq |f2) = j |f2) in units where h ^ I. 

In coarse-grained measurements our resolution is not able to resolve individual eigenvalues m of 
a spin component, say the z-component Jz, but bunches together Am neighboring outcomes into 
"slots" TO, where the measurement coarseness is much larger than the intrinsic quantum uncertainty 
of coherent states, i.e. Am^-^j. 

The question arises whether it is problematic to use coarse-grained (projective) von Neumann 
measurements of the form 



E (2-9) 



where |m) are the Jz eigenstates, as "classical measurements" as we did in the previous chapter. In 
contrast to the positive operator value measure (POVM), they have sharp edges and could violate 
the Leggett-Garg inequality by distinguishing with certainty between microstates at two sides of a 
slot border. E.g., if |to) and |to+1) belong to two different slots, macrorealism could be violated 
by a simple microscopic time evolution, producing cos{ujt) |to) + sm{ujt) |to+1). This is the reason 
why we model our coarse-grained Jz measurements as belonging to a (spin coherent state) POVM, 
where the element corresponding to the outcome to is represented by 



4tt 



fin. 



(2.10) 



Here, flm is the angular region of polar angular size AQm Am/j 3> 1/Vj whose projection onto 



the z-axis corresponds to the slot to (Figure 1.5 1. As the flm are mutually disjoint and form a 
partition of the whole angular region, we have 



^ — 



(2.11) 
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2 General conditions for quantum violation of macroscopic realism 



Due to overcompleteness of the spin coherent states the POVM elements are overlapping at the 
slot borders over the angular size which is small compared to the angular slot size A0m. 

To find out how Pfh looks like in the basis of eigenstates, we apply Pm to an arbitrary state 
|fc) and insert the identity operator 'Yl,k'W)(^'\ ~ ^- 



2j + l 



k' \ An 



{k'\n){n\k)d^n] \k' 



The scalar products in the above equation are of the form 

a/2 



cosJ+"f sin^'-™! e"™"^, 



(2.12) 



(2.13) 



where f2 = {•d^Lp). The azimuthal integration J^^^ exp[i(fc — A:')(^] d(p gives 2TrSkk' ■ The remaining 
integrand with fc = fc' is the modulus square |(fc|f2)p and independent of (/?, but we can just 
substitute the 27r factor from the performed i^s-integration by a trivial unperformed (/^-integration 
to keep the simple writing with an integration over flfn- Thus, 



2.1 + 1 



\{k\n)\'d'n \k){k\ 



(2.14) 



is diagonal in the basis of Jz eigenstates. 

The probability for getting the outcome ffi is given by 



Tr[pAr 



2j + l 
47r 



Inserting Pm from eq. (2.10), we obtain 



(2.15) 



2.1 + 1 
47r 



d^r^ / / d^r^' (n'\fi)(fi\p\fi') 

f2m 



(2.16) 



Using that JJd^Q \Q){n\ is the identity operator and renaming ft' to ft, we get 



J Jn^ 



This is exactly the integration of the positive probability distribution 

Q{n) ^ {n\p\n) 

associated to the quantum state p (the well-know Q-function [1]) over the region 



(2.17) 



(2.18) 



Given an arbitrary quantum state, under coarse-grained measurements the outcome 
probabilities at any time can exactly be computed from an ensemble of classical spins 
(i.e. there exists a hidden variable model). This is macrorealism per sc. 



2.3 Non-invasive measurability 
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2.3 Non-invasive measurability 

Upon a coarse-grained measurement with outcome m, the state p is reduced to 

. _ My-n p Mrn 
Prn — I 



(2.19) 



where we have chosen a particular (optimal) implementation of the POVM with the Hermitean 
Kraus operators 



Mm = Mi = J2 



/ 2j + l 
fe \l An 



\{k\n)\^d^n \k){k\ 



(2.20) 



It can be easily seen that they satisfy the necessary relation 



= P.l^ ■ 



We note that, independent of the concrete implementation, the Pm (and the Kraus operators) 
behave almost as projectors for all states except for those near a slot border: 



\Q) for Q. inside fim, 
for Q. outside fim. 



(2.21) 



In a proper classical limit {s/j /Am — s- 0) the relative weight of these border becomes vanishingly 
small. 

We now show that the Q-distribution before the measurement is the (weighted) mixture of the 
Q-distributions 

Q,^(r!) = {n\pm\n) (2.22) 

of the possible reduced states pm'- 



(2.23) 



This demonstrates that a fuzzy measurement can be understood classically as reducing 
the previous ignorance about predetermined properties of the spin system. 

The approximate sign reflects that, depending on the density matrix 



(2.24) 



this relationship may only approximately hold for the set of those CI = (?9, ip) near a slot border. 
In detail eq. ( |2.23| reads 



2j-H 
4 



E Cnn'{n\n){n'\n)^'-^J2 E Cnn'J2_V9Un)9rr^{n') {n\n){n'\n) (2.25) 



with 



gfh{k) 



2j + l 
47r 



\{k\n)Yd^n 



(2.26) 
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2 General conditions for quantum violation of macroscopic realism 




Figure 2.3: Left: The Q-distribution of the spin coherent state |f ,7r) located at the equator. Right: 
The mixture |(Q+ + Q~) of the Q-distributions of the two possible reduced states in a which-hemisphere 
measurement. Although this is a worst-case scenario in the sense that the quantum state is located exactly 
at the slot border, the two distributions are very similar. The pictures are drawn for the spin size j = 100 
but are scale invariant. 



which is smaller or equal to 1. The two sides are exactly equal if p is diagonal as J2fh 9fh{n) = Ifor 
all n. Deviations only occur if n, n' and j cos i? are all within a distance of order y/J to each other 
and to a slot border. If they are not close to each other, the quantity {U\n) {n' is exponentially 
small and suppression by the factor X]m\/ffni("') 9m{n') is not important. If n, n' are well within a 
slot, X]m\/5ni('^) 9m{n') IS almost identical to 1. 

Consider the worst-case scenario of a "which-hemisphere measurement" with only two slots 
m = + (northern hemisphere and fii — — (southern hemisphere f2_) and an initial spin 
coherent state 

Peq = |f,7r)(f,7r| (2.27) 

located exactly at the equator, i.e. the slot border. Both outcomes happen with the same probability 
Wj^ — w_ = i. Figure 2.3 illustrates the unmeasured state's Q-distribution, Q{^) = ^^i^ (ri|/5oq|ri). 



and the mixture \ [Q+{^) + Q_(ri)] of the Q-distributions 



Q^^^^^^l±l{m.P..M,\^) (2.28) 
Att 'w± 



of the reduced states, i.e. the left and right-hand side of eq. (2.23), respectively. Independent of 
the spin size j the Q-distributions with and without measurement have a very large overlap, 

j I \ [Q+(f]) + 0_(n)] Q(0) d^f] w 0.997 , (2.29) 



indicating the good quality of eq. (^2.23 1. Here, the overlap between two normalized probability 



distributions / and g is defined by // \/f{^) g{^) d^fi € [0, 1]. 



However, for non-invasiveness we need more than eq. (2.23). Consider the initial distribution 
of classical spins, Q(fl,tQ), corresponding to an initial quantum state p{to). We first compute the 
Q-distribution of the state p{tj) for an undisturbed evolution without measurement until some time 
tj, 

Q{n,t,)^'^^{n\p{t,)\n). (2.30) 



2.4 The sufficient condition for macroreahsm 
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Figure 2.4: Non-invasive measurability is fulfilled if the Q-distribution at some time tj, without any 
measurement before, is the same as the weighted mixture of the Q-distributions stemming from the possible 
reduced states from a measurement at some earlier time ti. 



This has to be compared with the mixture of all possible reduced distributions upon measurement 
at a time ti {to <ti < tj) with outcomes fh which evolved to tj, denoted as 

2j + 1 {n\Ut^.uMr-n P{ti) MrnUl-U I") 



(2.31) 



47r Wfn,U 

with Wfn^ti = Tr[/5(ti)Pm] the probability for outcome m at time ti and Ut = exp{—iHt) the time 



evolution operator (see Figure 2.4 1. The system evolves macrorealistically if these two quantities 
coincide for all ti and tj, i.e. if 



(2.32) 



This is the condition for non-invasive measurability (together with induction) . 



In a dichotomic scenario, the outcomes + and — correspond to finding the spin system in one out 
of two slots 771 = ±. This is represented by a measurement of two complementary regions ri+ and 
(for instance the northern and southern hemisphere in a "which hemisphere measurement"). 
Then, e.g., the probability for measuring — at t^ if + was measured at ti is given by 



93-11+ = // Q+,tA^,t3)'^^^ 
J Jf2_ 



(2.33) 



with Q-i-^tj (f2, is) the Q-distribution of the state which was reduced at ti with outcome -I- and 



evolved to ^3. If condition (2.32) is satisfied, it implies that the probabilities can be decomposed 
into "classical paths". This means that, e.g., 93_|i-)- is just the sum of the two possible paths via 
+ and — aX t2'. 

93-11+ = 92+I1+ 93-|2+,l+ + 92-I1+ 93-|2-,l+ , (2.34) 

where 93-|2+,i+ (93-|2-,i+) denotes the probability to measure — at ^3 given that -I- was measured 



at ti and + (— ) at t2- Thus, eq. (2.32 1 allows to derive Leggett-Garg inequalities such as ineq. (2.1 1 



2.4 The sufficient condition for macrorealism 

We can now establish the sufficient condition for macrorealism that holds even for isolated systems, 
namely 



An Ut m 



Ut \^) for one m, 

for all the others, 



(2.35) 
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2 General conditions for quantum violation of macroscopic realism 



for all t and fl, allowing deviations at slot borders. 



This means that the time evolution does not produce superpositions of macroscopically 
distinct states. 

Vice versa, if Ut produced states of the form a jfi') + (3\^") where neither \a\ nor is 
close to zero and with jfi') and jfi") belonging to macroscopically different outcomes (different 



slots), eq. (2.351 would not be fulfilled. Eq. (2.351 implies that P,n, and hence M^, quasi behave 

(2.36) 



as projectors and that 

{n\Ut^_u p{u) ul^u I") ~ E Wt^-uM^ p{u) Mr-nUl_t^ \n) 

This directly leads to eq. (2.32). Thus, eq. ( |2.35[ ) eq. (2.321 — ^ macrorealism. 



2.5 An oscillating Schrodinger cat 



We denote those Hamiltonians for which eq. (2.351 is satisfied under coarse-grained measurements 



as classical. An example is the rotation, say around x, Hiot — ^Jx, with the spin x-component 



and uj the angular precession frequency, which satisfies eq. (2.35) and moreover allows a Newtonian 



description of the time evolution as shown in the previous chapter. But there is no a priori reason 



why all Hamiltonians should satisfy eq. (2.35) 



Can one find non-classical Hamiltonians violating macrorealism despite coarse-grained 
measurements? 

The necessary condition for such a situation is that the Hamiltonian builds up coherences 
between states belonging to different slots. One explicit (extreme) example is 



H = iu{\-j){+j\-\+3){-3\) 



(2.37) 



which, given the special initial state |^(0)) = produces a time-dependent Schrodinger cat-like 
superposition of two distant (orthogonal) spin-j coherent states ('north') and ('south'): 



|*(t)) = cos(wi) + sm{Lut) 



(2.38) 



Under fuzzy measurements or apparatus decoherence after a premeasurement |104j . the state (2.38) 
appears like a statistical mixture at every instance of time: 



/5mix(i) = cos2(wi) \+j){+j\+sm^{ujt) . 



(2.39) 



While the two states Psup{t) = |^'(t))(^'(<)| and Pniix{t) have dramatically different P-functions and 
can be distinguished by sharp measurements, they are de facto equivalent on the coarse-grained 



level. The Q-distributions, Qsup for Psup{t) and Qmix for /5mix(^), are given by eq. (2.18). The 
coherence terms stemming from j5sup(t) are of the form (17 (— j| il) and vanish exponentially 
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fast with the spin length j for all fl. For j ^ 1 the Q-distributions are practically identical, i.e. 

Q — Qmix ~ Qsup- 

2j 



Q{n,t) = '^^-^ [cos2(c<jt)cos4^(^) +sin2(wt)cos4^(^)] , 



4tt 

where 01=1? (02 — tt ~ -d) is the angle between H. = {d, (p) and +z {—z). 
In general the P-function reads [TJ H] 

2j 



(2.40) 



^ PkqYkq{n,t) 

k=0 q=-k 



i-l)"-" ^{2j-k)\{2j + k+l)\ 



47r(2j)! 



Here, Ykq are the spherical harmonics and 

Pkq{t) = V2k + iy] i-iy 

' ^ m. 



^771, m—q 



it) 



3 k j 
-rn-\-q —qm 



(2.41) 



(2.42) 



where the last bracket denotes the Wigner 3j symbol. The c„„'(f) are the coefficients in the 



representation (2.24) of the density matrix one is interested in. 

The P and Q-functions of Psup and Pmix ^ = ^ are shown in Figure 2.5 for a certain choice of 



parameters. The P-function of the superposition is pathologically oscillating. The Q-distributions 
show just two peaks, corresponding to a classical mixture in which half of the spins are pointing 
into the north direction and the other into the south. Going to larger and larger values of j, 
i.e. from Schrodinger kittens to cats, makes it more and more difficult to observe the quantum 



nature of superposition states like (2.381. The angular resolution which is necessary to distinguish 



a superposition from the corresponding classical mixture is of the order of 1 /^/j. 



However, under the Hamiltonian (2.37 1 even simple dichotomic which-hemisphere measurements 



are sufficient to violate macrorealism. The temporal correlation function (for times ti and tj) reads 



Cij « cos[uj{tj-U)] 



(2.43) 



with an exponentially small correction due to the tiny chance that, e.g., \+j) can be found in 
the southern hemisphere. The system effectively behaves as a spin-^ particle which violates the 
Leggett-Garg inequality. It can be easily seen that without any measurement the initial state \+j) 
evolves to the state dX t2 — If, on the other hand, there is a measurement at ii = ^ 
in between, where we have an equal weight superposition of and the state is projected 

to either or |— j). Either of these reduced states will evolve to an equal weight superposition 
of and at ^2- Thus, eq. (2.321 is not fulfilled for fuzzy measurements and the evolution 



(2.371, as shown in Figure 2.6 



Despite the fact that apparatus decoherence or coarse-graining allow to describe the state 
effectively by a classical mixture at every instance of time, the non-classical Hamiltonian 
indeed builds up superpositions of macroscopically distinct states and allows to violate 
macrorealism. 



To get macrorealism one would have to coarse-grain always those states which are connected by 
the Hamiltonian (in time) and not necessarily in real space. In the present case it is (at least) the 
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Figure 2.5: Top left: The P-function Psup at time i = of the equal-weight superposition (2.38 1 of two 
opposite spin coherent states and |— j) for spin length j = 10, plotted in a rotated coordinate system 
in which = It is wildly oscillating with very large positive and negative regions. Top right: 



The P-function Pmix of the corresponding statistical mixture (2.391. Bottom: In every-day life the angular 
measurement resolution is much weaker than l/^/J (which is equivalent to Am^y^ in a measurement). 
Then we cannot distinguish anymore between the superposition state and the classical mixture, as both 
effectively lead to the same (positive) Q-distribution Qsup ~ Qmix- Nevertheless, the time evolution of such 
a mixture can violate macrorealism even under classical coarse-grained mesisurements. 



outcomes +j and — j that have to be coarse-grained into one and the same slot, which is of course 
highly counter-intuitive. Such a coarse-graining would lead to a different kind of macrorealistic 
physics than the classical laws we know, bringing systems through space and time continuously. 



2.6 Continuous monitoring by an environment 

Until now we have considered isolated systems. Under coarse-grained measurements or, mathemat- 
ically equivalent, apparatus decoherence (where the system is isolated and only after a premeasure- 
ment of the system an uncontrollable environment couples to the apparatus) , classical Hamiltonians 
lead to macrorealism whereas non-classical Hamiltonians allow to violate it. The question arises: 

What happens if the Hamiltonian is non-classical and the system is continuously moni- 
tored by an environment! 
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Figure 2.6: The time evolution of the Q-distribution for the non-classical (oscillating Schrodinger cat) 



Hamiltonian (2.371, shown for a spin j — 100. At time to — the quantum state is \+j) ('north') with the 
Q-distribution Q{to). At the later time ti ~ ^ the spin is in an equal-weight superposition of and 
and the Q-distribution Q{ti) shows two peaks at north and south. If no measurement takes place, 
the system reaches the state |— j) ('south') at t2 — represented by Q{t2). If, on the other hand, one 
performs a measurement at ti, the state is reduced to either or |— j), with Q-distributions Q+,ti{ti) 
and Q__tj(ii), respectively. Either of these two states will evolve into an equal-weight superposition until 
t2- The weighted mixture at that time, [Q±,ti (i2)+Q±,ti (t2)]/2 is different from the undisturbed evolution. 



Eqs. (2.351 and (2.32 1 are not fulfilled. Although at every instance of time the system can be described by 



a classical mixture, the time evolution of this mixture allows to violate the Leggett-Garg inequality and is 
in conflict with macrorealism. 



Given the non-classical Hamiltonian (2.371 with the time evolution operator 



Ut = +cos{L^t)i\+j)(+j\ + \-j){~j\) 

+ sm{c.t) (|-j)(+j|-|+j)(-j|) + 



(2.44) 



let us approximate the effects of system decoherence by the following simplified model: The initial 
state along north, 



P(0)-|+J>(+J|, 



(2.45) 
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freely evolves without decoherence a short time At to 

p{At) = m 

= cos'^{ujAt) + sin^(tjAi) \ + coh. terms , (2.46) 

where the coherence terms are of the form and Now we assume that the 

macroscopic spin system decoheres very rapidly (in the standard pointer basis of and |— j)), 
for instance due to the fact that a single qubit from the environment couples to it in a c-not 
manner |104j . becomes inaccessible immediately afterwards, and does not interact (recohere) with 
it anymore. If it is impossible to make (joint) measurements on the environmental qubit (and our 
spin system), the partial trace over the qubit of the total density matrix has to be performed, which 
kills the coherence terms in eq. (2.46), leading to the decohered state of the system: 

piAt) = cos2(c^At) + sin2(c^At) . (2.47) 

Assuming again free time evolution for a duration of At, this decohered state will evolve to 

p{2At) = UAt p{At) (2.48) 
= [cos^ {uj At) +sm^{uj At)] +2cos'^{ujAt)sm^{ijjAt) + coh. terms. 

Repeating the alternating sequence of rapid decoherence and free time evolution, we obtain the 
general expression for the (decohered) state at time nAt: 

p{nAt) = A,, + (1-A„) . (2.49) 

The survival probability to find the state along north. An, can be retrieved from the recurrence 
relation 

= 1 , (2.50) 

A„ = aA„_i + (l-a)(l-A.-i) (2.51) 

with integer n and 

a = cos^(cjAt) . (2.52) 

If At is not too small (to avoid a quantum Zeno-like freezing of the initial state |B6 ) but smaller 
than the dynamical timescale of the Hamiltonian, ^, the probability A„ decays to A^ = ^ in a 
way which can be very well approximated by 

A{t) = \{l-e-''') (2.53) 

with i the characteristic decay time. For t ^ ^ the state becomes an equal weight statistical 
mixture 

Pioo) . \+^)^+^\ + \-^)^-^\ (2.54) 

as illustrated in Figure |2.7[ 
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Figure 2.7: The survival probability An (black squares in red curve) in the decohered state p{nAt), 
eq. ( |2.49[ ), as a function of the number of steps, n. The free evolution time interval between instances 
of rapid decoherence is chosen to be At = j^^, where u is the angular frequency in the non-classical 
Hamiltonian ( 2.37 1. The state is driven into an equal weight mixture, eq. ( 2.54 1, with yloo — | (approaching 
the dashed line asymptotically). For comparison, we also draw the function cos^(tjt) — cos^ (ujnAt) which 
is the probability to find the state along north, i.e. \+j), at time t if no environmental decoherence takes 
place (dashed curve). 



The conclusion of this simple and crude model of decoherence — namely that the state is driven 
into a mixture where half of the spins point to north and half of the spins point to south — is 
expected to remain valid under more realistic circumstances where one does not separate into free 
evolution and rapid decoherence. As long as the environmental microscopic degrees of freedom only 
couple to the macroscopic spin system but do not disturb its diagonal elements, the system does 
not leave the subspace spanned by \+j) and and never populates any of the other states |m). 

Importantly, despite the non-classical Hamiltonian, the exponential decay of An due to deco- 
herence does not allow to violate the Leggett-Garg inequality any longer. If no (coarse-grained) 



measurement takes place, the spin's Q-distribution at time tj — i.e. the left-hand side of eq. (2.32 1 — 
is given by 

Q(tj) = A{t,) g„orth + [l-A{t,)] Qsouth, (2.55) 

where Qnorth (Qsouth) is the Q-distribution of a spin pointing to the north (south). If a measurement 
takes place at the intermediate time ti {0 < ti < tj)^ihe weighted mixture of the reduced and 



evolved Q-distributions — i.e. the right-hand side of eq. (^2.32 1 — takes on exactly the same form 



for all choices of ti and tj, given the exponential decay A{t) = ^(1— e~'^*), eq. (2.531, is used|^ 
Hence, the system's time evolution fulfills the condition ( |2.32 1 for non-invasive measurability, and 
consequently macrorealism is satisfied. 

However, decoherence cannot account for a continuous spatiotemporal description of the 
spin system in terms of classical laws of motion. 



^Note; After finding the spin along south at ti, which happens with probability 1—A{ti), A{tj —ti) is the (survival) 
probability to find the spin again along south at time tj . 

^It is interesting that any other form of the survival probability other than exponential decay violates the non- 



invasiveness condition (2.32 I 
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To see this, it is enough to use coarse-grained measurements corresponding to only three differ- 
ent angular regions, one covering the northern part, one the equatorial region, and one the southern 
part. The initial spin along north can be found pointing to the south at some later time, although 
it did not go through the equatorial region. No classical Hamilton function can achieve such discon- 
tinuous "jumps" of a spin vector. For claiming that classicality emerges from quantum physics, it 
is simply not enough to show that the density matrix is driven into a diagonal form due to tracing 
out the environmental degrees of freedom. One must also demonstrate that the time evolution of 
this mixed state can be described in terms of classical laws of motion. Classical physics is within 
the class of macrorealistic theories but it is more restrictive than macrorealism itself. 

The key point here is that in classical physics we have Newton's or Hamilton's differential equa- 
tions for observable quantities such as spin directions. Under all circumstances these equations 
evolve the observables continuously through real space. In quantum mechanics, however, the sit- 
uation is very different. The Schrodinger equation evolves the state vector continuously through 
Hilbert space but one cannot give a spatiotemporal description of the system's observables inde- 
pendent of observation. 

Figure |2.8| gives an overview of the conclusions we have reached until this point concerning 
the Leggett-Garg inequality, macrorealism, and classical laws of motion with respect to sharp and 
coarse-grained measurements as well as classical and non-classical Hamiltonians. 

2.7 Non-classical Hamiltonians are complex 

Finally, we suggest a possible reason why non-classical evolutions might be unlikely to be realized 
by nature: Such evolutions (i) either require Hamiltonians with many-particle interactions or (ii) 
a specific sequence of a large number of computational steps if only few-particle interactions are 
used. Then, they are of high computational complexity. In the first case the number of interacting 
particles, and in the second the number of computational steps has to scale linearly with the size 
of the Schrodinger cat state. Both cases intuitively seem to be of very low probability to happen 
spontaneously. 

Consider our spin-j as a macroscopic ensemble of iV = 2j spin-| particles (i.e. qubits) such 
as, e.g., any magnetic material is constituted by many individual microscopic spins. For violating 
macrorealism it is necessary to build up superpositions of two macroscopically distinct coherent 
states. For large j their angular separation A9 can be very small and only has to obey the coarse- 
graining condition A9 ^ l/VJ. This guarantees quasi-orthogonality as their (modulus square) 
overlap is cos {A0/2) - e-i'^'^'. Without loss of generality we consider again the particular 
Hamiltonian ( |2.37[ ). If |0) and |1) denote the individual qubit states 'up' and 'down' along z, then 
|ll...l) and |00...0) form the total spin coherent states and respectively. The Hamiltonian 
represents A'^-particle interactions of the form 

H = '-f {&®^ -af^) , (2.56) 

where a± = CT^i i o'y with ax and ay the Pauli operators. As an alternative one can simulate 
the evolution governed by this many-body interaction by means of a series of (in nature typically 
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Figure 2.8: Microscopic systems as well as macroscopic systems under sharp measurements allow to violate 
the Leggett-Garg inequality. Since no macroscopic "classical" observables are involved, one cannot speak 
about a violation of macrorealism in these cases. Under coarse-grained measurements of a macroscopic 
system and classical Hamiltonians not only macrorealism is valid but also classical laws of motion emerge. 
Non-classical Hamiltonians allow to violate macrorealism even under coarse-grained measurements. Deco- 
herence then establishes macrorealism but cannot account for a description of the system's time evolution 
in terms of classical laws of motion. 

appearing) few-qubit interactions (gates), using the methods of quantum computation science [67] . 
The task is to simulate 

I11...1) ^ cos(wt) I11...1) -f sin(w<) I00...0) . (2.57) 

Assuming next-neighbor qubit interactions, we start from the state |ll...l) and rotate the first qubit 
'1' by a small angle uiAt: |l)i — > cos(wAt) |l)i + sin(ijjAt) |0)i. Then we perform a controUed- 
not (c-not) gate between this qubit '1' and its neighbor '2' such that \x)i\ij)2 \x)i\xS)y)2 
{x,y — 0,1). Afterwards c-nots between qubits are performed such that all other qubits are 
reached. The whole procedure is depicted in Figure pT9| ^a) . This procedure brings us to the state 
at time At: — > cos(wAi) |ll...l) -I- sin(a;At) |00...0). To simulate the next time interval At, 
we have to undo all the c-nots, rotate the first qubit again by LuAt, and make all the c-nots again, 
leading to the correct state at time 2Ai. With this procedure we get a sequence of states, simulating 



the evolution (2.571. One needs 0{N) computational steps per interval At. This is known to be 



optimal in the case where only neighboring qubits can interact |15j . Relaxing this condition and 
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Figure 2.9: (a) In order to simulate the time evolution (2.571 of a qubit chain one has to rotate the first 



qubit by a small angle cuAt and sequentially make c-nots. For the next time interval At one has to undo the 
c-nots, rotate the first qubit again and make all the c-nots again. With this procedure one gets a sequence 



of states which approximate (2.571. (b) In contrast, the simulation of an interval At of a spin rotation of 
the whole chain can be achieved in a single global transformation on all qubits simultaneously. 



permitting two-qubit interactions between all possible qubits, allows to decrease the number of 
sequential steps but does not change the total number 0{N) of necessary gates per interval. 
Note for comparison, however, that the rotation (say around x) 



^rot=|V,' (2.58) 



2 



with k labeling the qubits, does not require multi-particle interactions. Moreover, the simulation 
of an interval At of a spin rotation of the whole chain, i.e. 

|111...) [cos(wAi) |1) +sin(tjAt) \0)f^, (2.59) 

can be achieved in a single global transformation on all qubits simultaneously, as shown in Fig- 



2.9 ». 



While both evolutions are rotations in Hilbert space (and require only polynomial re- 
sources), the simulation of the non-classical cosine-law between states that are distant 
in real space is — for macroscopically large N — computationally much more complex than 
the classical rotation in real space^ 



2.8 Information and randomness 

In a perfect von Neumann measurement of a spin component with sharp resolution an individual 
state I to) out of the 2j + 1 w 2j possible ones carries 

/sharp « l0g2(2j) = 1 + \og^ ] (2.60) 

bits of information. Coarse-grained measurements correspond to the fact that we cannot resolve 
individual eigenvalues to but only whole bunches of size Ato 3> \/~]. The finding that an outcome 
lies in a certain slot of size Ato = c\fj (with c 3> 1) carries only 

/e.g. « log2(^) = 1 - log2 C + \ log2 J (2.61) 

^Onc should, however, mention the possibility that an external field may produce an effectively simple non-classical 
Hamiltonian for N qubits where the field interacts with the collective modes = |ll...l) and = |00...0). 
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bits of information. For large j, i.e. j 3> c 3> 1, the information gain in a sharp quantum mea- 
surement is approximately logj j bits, whereas in the classical case it is (at most) only half of that, 
namely ^ logj j bits [5S] : 

^c.-g. ~ 2 -^sharp ■ (2.62) 

Finally, we note that — given coarse-grained measurements — it is objectively random which of the 
two states, 'north' or 'south', one will find in a spin measurement in the Schrodinger cat state (2.38 1. 
Classical physics emerges out of the quantum world but the randomness in the classical mixture 
is still irreducible. Which possibility becomes factual is objectively random and does not have a 
causal reason, according to the Copenhagen interpretation. 

In the last chapter we will address the question of quantum randomness in more detail and link 
it with mathematical undecidability. 
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Chapter 3 



Entanglement between 
macroscopic observables 

Summary: 

Wc investigate entanglement between collective operators of two blocks of oscillators in an infi- 
nite linear harmonic chain. These operators are defined as averages over local operators (individual 
oscillators) in the blocks. On the one hand, this approach of "physical blocks" meets realistic 
experimental conditions, where measurement apparatuses do not interact with single oscillators 
but rather with a whole bunch of them, i.e. where in contrast to usually studied "mathematical 
blocks" not every possible measurement on them is allowed. On the other, this formalism natu- 
rally allows the generalization to blocks which may consist of several non-contiguous regions. We 
quantify entanglement between the collective operators by a measure based on the Peres-Horodecki 
criterion and show how it can be extracted and transferred to two qubits. Entanglement between 
two blocks is found even in the case where none of the oscillators from one block is entangled with 
an oscillator from the other, showing genuine bipartite entanglement between collective operators. 
Allowing the blocks to consist of a periodic sequence of subblocks, we verify that entanglement 
scales at most with the total boundary region. We also apply the approach of collective operators 
to scalar quantum field theory. 

What can we learn about entanglement between individual particles in macroscopic samples by 
observing only the collective properties of the ensembles? Using only a few experimentally feasible 
collective properties, we establish an entanglement measure between two samples of spin-i particles 
(as representatives of two-dimensional quantum systems). This is a tight lower bound for the 
average entanglement between all pairs of spins in general and is equal to the average entanglement 
for a certain class of systems. We compute the entanglement measures for explicit examples and 
show how to generalize the method to more than two samples and multi-partite entanglement. On 
the fundamental side, our method demonstrates that there is no reason in principle why purely 
quantum correlations could not have an effect on the global properties of objects. On the practical 
side, it enables us to characterize the structure of entanglement in large spin systems by performing 
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only a few feasible measurements of their collective properties, independently of the symmetry and 
mixedness of the state. 

Since this analysis uses sharp measurements, it is not in disagreement with our quantum-to- 
classical approach resting upon coarse-grained measurements. 

This chapter mainly bases on and also uses parts of References [531 HI] • 

• J. Kofler, V. Vedral, M. S. Kim, and C. Brukner 

Entanglement between collective operators in a linear harmonic chain 
Phys. Rev. A 73, 052107 (2006). 

• J. Kofler and C. Brukner 

Entanglement distribution revealed by macroscopic observations 
Phys. Rev. A 74, 050304(R) (2006). 
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3.1 The linear harmonic chain 

Quantum entanglement is a physical phenomenon in which the quantum states of two or more 
systems can only be described with reference to each other, even though the individual systems 
may be spatially separated. This leads to correlations between observables of the systems that 
cannot be understood on the basis of classical (local realistic) theories [5]. Its importance today 
exceeds the realm of the foundations of quantum physics and entanglement has become an important 
physical resource, like energy, that allows performing communication and computation tasks with 
efficiency which is not achievable classically [67 . Moving to higher-dimensional entangled systems 
or entangling more systems with each other, will eventually push the realm of quantum physics well 
into the macroscopic world. It will therefore be important to investigate under which conditions 
entanglement within or between "macroscopic" objects, each consisting of a sample containing a 
large number of the constituents, can arise. 

Recently, it was shown that macroscopic entanglement can arise "naturally" between con- 
stituents of various complex physical systems. Examples of such systems are chains of interacting 
spin systems [H71 |5j , harmonic oscillators [3 [SS] and quantum fields [5T] . Entanglement can have an 
effect on the macroscopic properties of these systems [371 ISHl [IH] and can be in principle extractable 
from them for quantum information processing [5TJ [7TJ HOI 113 • 

With the aim of better understanding macroscopical entanglement we will investigate entangle- 
ment between collective operators. A simple and natural system is the ground state of a linear chain 
of harmonic oscillators furnished with harmonic nearest-neighbor interaction. The mathematical 
entanglement properties of this system were extensively investigated in [TJ |T3J ESI [H] • Entangle- 
ment was computed in the form of logarithmic negativity for general bisections of the chain and 
for contiguous blocks of oscillators that do not comprise the whole chain. It was shown that the 
log-negativity typically decreases exponentially with the separation of the groups and that the 
larger the groups, the larger the maximal separation for which the log-negativity is non-zero [7]. 
It also was proven that an area law holds for harmonic lattice systems, stating that the amount of 
entanglement between two complementary regions scales with their boundary |26l llOOj . 

In a real experimental situation, however, we are typically not able to determine the complete 
mathematical amount of entanglement (as measured, e.g., by log-negativity) which is non-zero even 
if two blocks share only one arbitrarily weak entangled pair of oscillators. Our measurement appa- 
ratuses normally cannot resolve single oscillators, but rather interact with a whole bunch of them 
in one way, potentially even in non- contiguous regions, thus measuring certain global properties. 
Here we will study entanglement between "physical blocks" of harmonic oscillators — existing only 
if there is entanglement between the collective operators defined on the entire blocks — as a function 
of their size, relative distance and the coupling strength. Our aim is to quantify (experimentally 
accessible) entanglement between global properties of two groups of harmonic oscillators. Surpris- 
ingly, we will see that such collective entanglement can be demonstrated even in the case where 
none of the oscillators from one block is entangled with an oscillator from the other block (i.e. it 
cannot be understood as a cumulative effect of entanglement between pairs of oscillators), which 
is in agreement with Reference j^. This shows the existence of bipartite entanglement between 
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Figure 3.1: In a linear iiarmonic cliain tiie oscillators are situated in a harmonic potential with frequency 
uj. Each oscillator, labeled with the index n, is coupled with its nearest neighbors, n ~ 1 and n + 1, by a 
harmonic potential with the coupling frequency fl. 



collective operators. 

Because of the area law llOOj the amount of entanglement is relatively small in the first 
instance. We suggest a way to overcome this problem by allowing the collective blocks to consist 
of a periodic sequence of subblocks. Then the total boundary region between them is increased 
and we verify that indeed a larger amount of entanglement is found for periodic blocks, where the 
entanglement scales at most with the total boundary region. We give an analytical approximation 
of this amount of entanglement and motivate how it can in principle be extracted from the chain |81l 

iniiHniEi!. 

Methodologically, we will quantify the entanglement between collective operators of two blocks 
of harmonic oscillators by using a measure for continuous variable systems based on the Peres- 
Horodecki criterion [75J HU |HH1 The collective operators will be defined as sums over local 
operators for all single oscillators belonging to the block. The infinite harmonic chain is assumed 
to be in the ground state and since the blocks do not comprise the whole chain, they are in a mixed 
state. 

We investigate a linear harmonic chain, where each of the N oscillators is situated in a harmonic 
potential with frequency to and each oscillator is coupled with its neighbors by a harmonic potential 



with the coupling frequency f2 (Figure 3.11. The oscillators have mass m and their positions 
and momenta are denoted as and p^, respectively. Assuming periodic boundary conditions 
Oin+i = the Hamilton function reads |84j 




2 7^2 ^ r>2 (-77 77 ^2 \ 

(3.1) 



mu)^ qj {qj — qj_i) 

2 ^ 2 




We canonically introduce dimensionless variables: qj = C q^ and pj = Pj/C, where C is given by 
C = ^/mu;il + 2n^/u;-^y/^ [H. By this means the Hamilton function becomes 



(3.2) 



with the abbreviations a = 2 51^/(217^ +uj'^) and Eq = y/2 il'^ + uj'^. The (single) coupling constant 
is restricted to values < a < 1, where a — > in the weak coupling limit (f2/w — > 0) and a — > 1 in 
the strong coupling limit (fl/uj oo). 

In the language of second quantization the positions and momenta are converted into operators 
(qj qj, Pj — > Pj) and are expanded into modes of their annihilation and creation operators, a 
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Figure 3.2: The two-point vacuum correlation functions (3.6 1 and (3.71 for the linear harmonic chain with 
coupling a = 0.5. The blue curve shows = {0| Qi qj 

with I = \i — j\- Only integer values for / are meaningful. 



0), and the red curve shows h^i. 
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Here dk = li^kjN (with fc = 0, 1, ...,N — 1) is the dimensionless pseudo-momentum and 

v(Qu) = V 1 - " cos 



(3.3) 
(3.4) 

(3.5) 



is the dispersion relation. The annihilation and creation operators fulfil the well known commuta- 
tion relation [a((?fe), a^(^?fc')] = ^kk' , since = i(5ij has to be guaranteed. The ground state 
(vacuum), denoted as |0), is defined by aid]^ |0) = holding for all 6*^,. The two-point vacuum 
correlation functions are 



*1»-Jl 



|0) = {v^h) - (2^)" E,^o ^(^fc)cos(Z^fc), 



(3.6) 
(3.7) 



where I = \i — i\. In the limit of an infinite chain (A^ — s- oo) — which we will study below — 
and for I < N/2 they can be expressed in terms of the hypergeometric function 2F1 [E]: gi = 
[7J/{2fi)]{'-y^) 2Fi(l/2, 1 + 1/2, 1 + 1, z2), hi = (^z72) {'-''/^) 2i^i(-l/2, 1 - 1/2, 1 + 1, z^), where 
z = (1 - Vl - a^)/a and = l/VIT^ (Figure Eil. 



3.1.1 Defining collective operators 

In the following, we are interested in entanglement between two "physical blocks" of oscillators, 
where the blocks are represented by a specific form of collective operators which are normalized sums 
of individual operators. By means of such a formalism we seek to fulfil experimental conditions 
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B 



Figure 3.3: Two blocks A and B of a harmonic chain. Each block consists of n oscillators and the blocks 
are separated by d oscillators. 



and constraints, since finite experimental resolution implies naturally the measurement of, e.g., 
the average momentum of a bunch of oscillators rather than the momentum of only one. On the 
other hand, this formahsm can easily take account of blocks that consist of non- contiguous regions, 
leading to interesting results which will be shown below. We want to point out that this convention 
of the term block is not the same as it is normally used in the literature. In contrast to the latter, 
for which one allows any possible measurement, our simulation of realizable experiments already 
lacks some information due to the averaging. 

Let us now consider two non-overlapping blocks of oscillators, A and B, within the closed 
harmonic chain in its ground state, where each block contains n oscillators. The blocks are separated 



by 0? > oscillators (Figure 3.3 1. We assume n,d ^ N and iV — > cx) for the numerical calculations 
of the two-point correlation functions. 

By a Fourier transform we map the n oscillators of each block onto n ("orthogonal") frequency- 
dependent collective operators 

I 27rijfc 



9, e n , (3.8) 



,(fc) _ 1 ^r^ . -^TTijfc 



P('=)^^^p,e-— , (3.9) 
V jeA 

with the frequencies fc = 0, n — 1, and analogously for block B. The commutator of the collective 
position and momentum operators is 

[Qf,Pp]^iSkk'. (3.10) 

This means that collective operators for different frequencies k ^ k' commute. For different blocks 
the commutator vanishes: [<5a ^ -^s^ ~ ^■ 

If the individual positions and momenta of all oscillators are written into a vector 

X = {qi,pi,q2,P2, ...,qN,PN)^, (3.11) 

with T denoting the transpose, then there holds the commutation relation 

[xi,Xj]=inij (3.12) 

with r2 the n-fold direct sum of 2 x 2 symplectic matrices: 

(3.13) 
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A matrix S transforms x into a vector of collective (and uninvolved individual) oscillators: 

±^Sic=i{Qf,Pi'^,,{Q^B\Pi^}k,{qs,Pshf- (3.14) 

Here {(5i^\ Pi'^^j/c = (Q^°\ P4°\ Qa^^\Pa^^^) denotes all collective oscillators of block A and 
analogously for block B, whereas {qj,Pj}j denotes the 2 (N — 2n) position and momentum entries 
of those N — 2n oscillators which are not part of one of the two blocks. The matrix S corresponds 
to a Gaussian operation [31' . It has determinant det S = 1 and preserves the symplectic structure 

rj = S^rJS, (3.15) 

and hence 

[X,,X^]^in,j (3.16) 



for all in particular verifying eq. (3.101. This means that the Gaussianness of the ground 
state of the harmonic chain (i.e. the fact that the state is completely characterized by its first and 
second moments) is preserved by the (Fourier) transformation to the frequency-dependent collective 
operators. 

3.1.2 Quantifying entanglement between collective operators 

In reality, we are typically not capable of single particle resolution measurements but only of 
measuring the collective operators with one frequency, namely k = 0, i.e. the "average" over the 
individual oscillators. Note that in general the correlations of higher-frequency collective operators, 
e.g., HQa^)'^) or (Q^^Q^"*) with fc 7^ 0, are not real numbers. Therefore, as a natural choice, we 
denote as the collective operators 



{A 

1 / n 

j<^A 



1 (3.17) 



/n 
jeA 



and analogously for block B. It seems to be a very natural situation that the experimenter only has 
access to these collective properties and we are interested in the amount of (physical) entanglement 
one can extract from the system if only the collective observables Qa,b and Pa,b are measured. 

Reference [88] derives a separability criterion which is based on the Peres-Horodecki criterion [75j 
|44] and the fact that — in the continuous variables case — the partial transposition allows a geometric 
interpretation as mirror reflection in phase space. Following largely the notation in the original 
paper, we introduce the vector 

e= (Qa,A4,Qs,Ps) (3.18) 

of collective operators. The commutation relations have the compact form [^q,^/?] = iKap with 
K = ®^^]^( J^i ) . The separability criterion bases on the covariance matrix (of first and second 



7j = H^-l 0^ 

moments) 

Vc.p=^{ A^M^ + A^ACa ) , (3.19) 
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where A^a = ~ (Ca) with (^q,) = in our case (the ground state is a Gaussian with its mean at 
the origin of phase space). 

The covariance matrix V is real (which would not be the case for higher-frequency collective 
operators) and symmetric: (QaQb) = {QbQa) and {PaPb) = {PrPa), coming from the fact 



that the two-point correlation functions (3.6 1 and (3.7 1 only depend on the absolute value of the 
position index difference. On the other hand, using eqs. ( |3.3[ ) and (3.4|, we verify that ((jiPj) = 
i(2^)"'Ef=o'exp[i0fe(i - j)] and {p, q,) = -i (2 7V)-i Efjo' exp[i 0fe(j - t)]. For i ^ j both 
summations vanish {0^. = 2nk/N and i,j integer) and for i = j they are the same but with 



opposite sign. Thus, in all cases {qiPj\ 
operators and hence we obtain 



{pj Qi). These symmetries also hold for the collective 



(3.20) 



The matrix elements are 

(3.21) 

_ (3.22) 

-_A 

jeA ieB 

= ^EE%-I- (3-24) 

jeA ieB 

To quantify entanglement between two collective blocks we use the degree of entanglement e, given 
by the absolute sum of the negative eigenvalues of the partially transposed density operator: e = 
Tr|pT«|-l, i.e. by measuring how much the mirror reflected state p^^ , where the momenta in block 
B are reversed, fails to be positive definite. This measure (proportional to the negativity) is based 
on the Peres-Horodecki criterion [751 31] and was shown to be an entanglement monotone [5^ I93j . 
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[PaPb] 



For covariance matrices of the form (3.201 it reads [50 



e = max | 0, — — 1 

O1O2 



(3.25) 



with 

6i = G~\Gab\, 
S2 = H-\Hab\- 

In general, the numerator is defined by the square of the Heisenberg uncertainty relation 

{6162)0= \{[Qa.b,Pa,b])\'^\, 



(3.26) 
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d=0 d=l 




Figure 3.4: Degree of collective entanglement e for two blocks of oscillators as a function of their size 
n. Left: The blocks are neighboring (d — 0) and entanglement exists for all n and coupling strengths a. 
Plotted are a — 0.99 (red), a — 0.9 (green) and a = 0.5 (blue). Right: The same for two blocks which are 
separated by one oscillator (d = 1). The two blocks are unentangled for n = 1 but can be entangled, if one 
increases the block size (n > 1), although none of the individual pairs between the blocks is entangled. 



where the last equal sign holds due to eq. (3.101. We note that e is a degree of entanglement (in 



the sense of necessity and sufficiency) only for Gaussian states which are completely characterized 
by their first and second moments, as for example the ground state of the harmonic chain we are 
studying. However, we left out the higher-frequency collective operators (and all the oscillators 
which are not part of the blocks) and therefore, the entanglement e has to be understood as the 
Gaussian part of the amount of entanglement which exists between (and can be extracted from) 



the two blocks when only the collective properties Qa,b and Pa,b, as defined in eq. (3.171, are 
accessible. 

There also exists an entanglement witness in form of a separability criterion based on variances, 
where 

^ = {{Qa - Qb?) + {{Pa + Pb?) = "^{0 - Gab + H + Hab) <2 (3.27) 
is a sufficient condition for the state to be entangled We note that the above negativity 



measure (3.251 is "stronger" than this witness in the whole parameter space (a,n). In particular, 
there are cases where e > although A > 2. This is in agreement with the finding that the variance 
criterion is weaker than a generalized negativity criterion |86j . 



We further note that the amount of entanglement ( |3.25 1 is invariant under a change of potential 



redefinitions of the collective operators, e.g., Qa = '^j^A'lj "-"^ '9-4 = ( V^*) EjeA '?j ' ^ then the 
modified scaling in the correlations (G, Gab, H, and Hab) is exactly compensated by the modified 
scaling of the Heisenberg uncertainty in the numerator. 



Figure 3.4 shows the results for d = and d=l. In the first case — if the blocks are neighboring — 
there exists entanglement for all possible coupling strengths a and block sizes n. In the latter case — 
if there is one oscillator between the blocks — due to the strongly decaying correlation functions g 
an h there is no entanglement between two single oscillators (n = 1), but there exists entanglement 
for larger blocks (up to n = 4, depending on a). The statement that entanglement can emerge by 
going to larger blocks was also found in [7]. But there the blocks were abstract objects, containing 
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Figure 3.5: Two periodic blocks of a harmonic chain A and B. Each block can consist of m subblocks 
with s oscillators each, separated by d oscillators. In the picture d=l, m = 2, s = 3 and the number of 
oscillators per block is n = ms — 6. 



all the information of their constituents. In the case of collective operators, however, increasing 
the block size (averaging over more oscillators) is also connected with a loss of information. In 
spite of this loss and the mixedness of the state, two blocks can be entangled, although none of the 
individual pairs between the blocks is entangled — indicating true bipartite entanglement between 
collective operators (multipartite entanglement between individual oscillators). For d > 2, however, 
no entanglement can be found anymore. 

These results are in agreement with the general statement that entanglement between a region 
and its complement scales with the size of the boundary [26l llOOj . In the present case of two 
blocks in a one-dimensional chain (Figure 3.3 1 the boundary is constant and as the blocks are 
made larger, the entanglement decreases since it is distributed over more and more oscillators. We 
therefore propose to increase the number of boundaries by considering two non-overlapping blocks, 
where we allow a periodic continuation of the situation above, i.e. a sequence of to > 1 subblocks. 



separated by d oscillators and each consisting of s > 1 oscillators, where ms — n (Figure 3.5). 

The degree of entanglement between two periodic blocks of non-separated (d = 0) one-particle 
subblocks (s = 1) is larger for stronger coupling constant a and grows with the overall number of 
oscillators n (Figure 3.6 1). For given a and n and no separation between the subblocks {d — 0) the 
entanglement is larger for the case of small subblocks, as then there are many of them, causing a 



large total boundary (Figure 3.6 d). Entanglement can be even found for larger separation {d = 1, 2) 
with a more complicated dependence on the size s of the subblocks. There is a trade-off between 
having a large number of boundaries and the fact that one should have large subblocks as individual 
separated oscillators are not entangled (Figure 3.6 :,d). For c? > 3 no entanglement can be found 
anymore. (In a realistic experimental situation, where the separation d is not sharply defined, e.g., 
where there are weighted contributions for d = 0, 1, dmax, entanglement can persist even for 
dmax > 3, depending on the weighting factors.) 

For the sake of completeness we give a rough approximation of the entanglement between two 
periodic blocks. Let us assume that the subblocks are directly neighbored, d ^ 0. Furthermore, we 
consider couplings a such that we may neglect higher than next neighbor correlations (a < 0.5), 
i.e. we only take into account go, gi, ho and hi. The correlations read 

^ = - E E ^90+ \ 91 , (3.28) 

^ 1 \ ^ \ - 2to — 1 ,„ 

GAB = -Y.T.9U-r\^^^9i, (3.29) 
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Figure 3.6: Degree of collective entanglement e for two periodic bloclcs of oscillators as a function of their 
total size n. (a) Neighboring one-particle subblocks {d = 0, s — 1). Entanglement monotonically increases 
with n and becomes larger as the coupling strength a increases. Plotted are a = 0.99 (red), a = 0.9 (green) 
and a = 0.5 (blue), (b) The coupling is fixed to a = 0.99 for this and the subsequent graphs. There is 
no separation, d = 0. Plotted are the cases s = 1,2,5 in red, green, blue, respectively. For fixed n the 
entanglement is more or less proportional to the number of boundaries, i.e. inversely proportional to the 
subblock size s. (c) and (d) correspond to the cases d — 1 and d = 2, respectively. The dependence on the 
size of the subblocks is more complicated as there is a trade-off between having a large number of boundaries 
(i.e. small s) and the fact that one should have large subblocks as individual separated oscillators are not 
entangled. 



and analogously for H and Hab- The first equation reflects that there are n self-correlations and 
m{s — 1) nearest neighbor pairs (which are counted twice) within one block, i.e. s — 1 pairs per 
subblock. The second equation represents the fact that there are 2m — 1 boundaries where blocks 



A and B meet. Using s = n/m, the entanglement (3.251 becomes (note that gi > and hi < 0) 

1 



- 1. 



(3.30) 



4[5o + (2-^)5i] N + (2-i)/Ji] 

For given n this approximation obviously increases with the total number of boundaries, m. It can 
be considered as an estimate for a situation like in Figure 3.6 3, if a smaller coupling is used such 
that the neglect of higher correlations becomes justified. 



We close this section by annotating that the entanglement (3.251 between collective blocks of 
oscillators — being the Gaussian part — can in principle (for sufficient control of the block separation 
d) be transferred to two remote qubits via a Jaynes-Cummings type interaction [5TJ |7TJ |5D] . For 
the interaction with periodic blocks "gratings" have to be employed in the experimental setup. The 
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interaction Hamiltonian is of the form 



2 



(3.31) 



where uii is the Rabi frequency and &^ = {a- ^ = \e)ii{g\ is the bosonic operator (with \g)i and 
\e)i the ground and the excited state) of the i-th qubit {i = 1,2). 



3.1.3 Collective operators for scalar quantum fields 

The continuum limit of the linear harmonic chain is the (l+l)-dimensional Klein-Gordon field 
(j){x^t) with the canonical momentum field ■n{x,t) — (j)(x,t). It satisfies the Klein-Gordon equation 
(in natural units h = c = I) with mass m 



(3.32) 



With the canonical quantization procedure ((> and tt become operators satisfying the non-trivial 
commutation relation [(/)(x, i), 7r(x', i)] = id{x — x'). The field operator can be expanded into a 
Fourier integral over elementary plane wave solutions |12j 

dfc 



(j){x,t) = 
7r(a;, t) = — i 



dfcwfc 



a{k) e''=^-'"'=* H.c] , 
[a(/c) e"=^-'"''* - H.c] , 



(3.33) 
(3.34) 



where k is the wave number and ojk — +\/k'^ + m? is the dispersion relation. The annihilation and 
creation operators fulfil [a(A;), a^(fc')] — 5{k — k'). We write the field operator as a sum of two con- 
tributions <b = <h^^'> + 



where 



(-) 



) is the contribution with positive (negative) frequency. 



Thus, (^'+^ corresponds to the term with the annihilation operator in eq. (3.331. The vacuum corre- 
lation function is given by the (equal-time) commutator of the positive and the negative frequency 
part: 

(O|0(x,t)0(y,t)|O) = [^(+)(a;,t),^(-)(y,t)]. (3.35) 

It is a peculiarity of the idealization of quantum field theory that for x = y this propagator diverges 
in the ground state: 



\i\x,t)\Q) 



(3.36) 



The same is true for (0| 7r^(x,t) |0) and hence we cannot easily build an entanglement measure like 



for the harmonic chain, since the analogs of the two-point correlation functions go and /iq, eqs. (3.6) 



and (3.7 1, are divergent now. Automatically, we are motivated to study the more physical situation 
and consider extended space-time regions, which means that we should integrate the field (and 
conjugate momentum) over some spatial area. We define the collective field operators 







fL/2 


h Xq, t) dx , 






'-L/2 






rL/2 








/ Tt{x - 
'-L/2 


- Xq, t) dx . 



(3.37) 
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Therefore, ^L(xo,t) and IlL(xo,t) are equal-time operators which are spatially averaged over a 
length L, centered at position a:o- The commutator is 



1 /-^/^ /-^/^ 

[^L{xa,t),IlL{xo,t)] = - / iS{x-y)dxdy = i, 

^ J-L/2 J-L/2 



(3.38) 



which is in complete analogy to eq. (3.101. If <i>L and 11^ correspond to separated regions without 



overlap, i.e. |xo — 2/o| > L, then of course [<&l(xo, t), nL(2/0; i)] — 0. The spatial integration in 

sin(^) [a{k) e"=^o-''^'=* + H.c] , 



eq. (3.371 can be carried out analytically: 

dfc 



1 



/ttL 

— i 



dk ^Jujk 



sm 



(^) [a(fc)e"=^°-'"''*-H.c.]. 



(3.39) 
(3.40) 



V TtL j-oo 

The final step is to calculate the propagators of the field and the conjugate momentum. We find 



\\<i>L{xo,t)<i>L{yo,t)\0) 
1 f°° dk 

= / : SlI 

nL k'^y'k'^ + vn? 
i^n_^(r) = (0|nL(xo,t)ni(yo,i)|0) 
1 dk^/}? ^ 



2fkL 



) cos(fc r) 



,„ — sin^(4r)cos(fcr) , 



(3.41) 



(3.42) 



with r = \xq — yol the distance between the centers of the two regions, reflecting the spatial sym- 
metry. Thus ^(0) and Df^ ^(0) are the analogs of {Q\ g) and {P\ g) (mfra-block correlations 
within the same block), respectively, whereas j^{r > L) and Df^ j^{r > L) correspond to ( QaQb ) 
and {PaPb) (mier-block correlations between separated blocks). 

The expressions ( |3.41 1 and (3.421 are finite, especially for r = 0. Mathematically, the integration 
over a finite spatial region L corresponds to a cutoff, which removes the divergence we faced in 

0. 



eq. (3.361. However, the expressions are ill defined for L 



Applying the entanglement measure (3.25) with G 



(0), Gab = D^^r), 

and Hab = Djj ^(r) does not indicate entanglement for any choice of L and r > L. The same is 
true for the generalized case of blocks consisting of periodic subregions of space, showing an inherent 
difference between the harmonic chain and its continuum limit. This might be due to the fact, that 
any spatial integration immediately corresponds to an infinitely large block in the discrete harmonic 
chain and that the information loss (compared to the mathematical indeed existing exponentially 
small entanglement [gT]) due to the collective operators already is too large. Nonetheless, defining 



collective operators like in eq. (3.371 and use of the measure (3.251 may reveal entanglement between 



spatially separated regions for other quantum field states. 



3.2 Spin ensembles 

Observation of quantum entanglement between increasingly larger objects is one of the most promis- 
ing avenues of experimental quantum physics. Eventually, all these developments might lead to a 
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full understanding of the simultaneous coexistence of a macroscopic classical world and an under- 
lying quantum realm. Macroscopic samples typically contain N ~ 10^° particles. Because the 
system's Hilbert space grows exponentially with the number of constituent particles, a complete 
microscopic picture of entanglement in large systems seems to be in general intractable. The ques- 
tion arises: What can we learn about entanglement between constituent particles of macroscopic 
samples, if only limited experimentally accessible knowledge about the samples is available? 

There is a strong motivation in addressing this question because of recent experimental progress 
in creating and manipulating entangled states of increasing complexity, such as spin-squeezed states 
of two atomic ensembles |47j . In such experiments one typically measures only expectation val- 
ues of collective operators of two separated samples. It is known that the two samples of spins 
can be characterized as either entangled or separable by measuring collective spin operators 89J. 
Furthermore, such measurements are shown to be sufficient to determine entanglement measures 
of Gaussian states [ST] and of a pair of particles that is extracted from a totally symmetric spin 
state (invariant under exchange of particles) [SB]- It appears that collective operators cannot be 
used to fully characterize entanglement in composite systems without strong requirements on the 
symmetry of the state. 

Here we present a general and practical method for entanglement detection between two samples 
of spins. It solely employs collective spin properties of the samples and works irrespectively of the 
number of spin particles constituting the samples and with no assumption about the symmetry 
or mixedness of the state. The method is based on an entanglement measure which is a tight 
lower bound for average entanglement between all pairs of spins belonging to the two samples. This 
measure is equal to the average entanglement for a certain class of systems which need not be totally 
symmetric. We generalize the method to obtain the entanglement measure between M separated 
spin samples based on collective measurements. The results apply for any entanglement monotone 
that is a convex measure on the set of density matrices (e.g., concurrence [H], negativity |105ll93j . 
three-way tangle [21] )• 



3.2.1 Bipartite entanglement 



Consider an ensemble of spin- ^ particles which is separated into two ensembles A and B (Figure 3.7 1. 
Each of these ensembles contains a large number of spins which we denote as n. Because of the 
large dimensions, 

d = 2" , (3.43) 

of the samples' Hilbert spaces, the structure of entanglement between the two samples is consider- 
ably more complex than between two single spins. While there are experimentally viable methods 
for detecting entanglement, they still require a large number of parameters to be determined (pro- 
portional to (P [15]V 

The problem simplifies in situations in which each of the ensembles of n spins can be treated 
as one large total spin of length ^. This means that, within the ensembles, the individual spin-^ 
particles form symmetrized states (Dicke states) . Though this reduces the dimension of the Hilbert 
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A B 




Figure 3.7: Two contiguous and non-overlapping spin subsystems A and B, each of which contains a large 
number n of spins. What can we learn about entanglement of a pair (a, /3) of spins chosen at random where 
a £ A and /3 € -B, if individual spins are experimentally not accessible but only the collective properties 
of the samples A and B? What can we learn about entanglement between A and B from such collective 
measurements? 



space of A (B) to 



d' 



1, 



(3.44) 



entanglement determination is still demanding for large n both experimentally and theoretically: 
Analytical solutions exist only for pure states in general and for mixed states only for small n 82J . 

In this section, we give a method to detect entanglement between large spin samples by mea- 
suring only a small number of collective spin properties (sample spin components and their corre- 
lations), which is independent of the sample size n. The collective spin operators are 



i I 



(/3) 



(3.45) 



The index i denotes the spatial component of the spins: i G {l = x,2 = y,3 = z}. The Pauli matrix 
of the spin at site a e A is given by (t|"\ and analogously for a spin (3 from subsystem B. Note 
that the collective operators satisfy the usual commutation relations 



(3.46) 



smce [al \aj '\ 



2i5, 



The spin expectation values and correlations are 



and analogously for Sf ~ ^^peB 



2 fc2 2 fc2 , 

n h n h ] 



5.(/3), where 5/^,5f e[-f ,f], 7;j^e[ , , , 



(3.47) 
(3.48) 

These 



are only 15 numbers. Here gi{a), gi{(3) and hij(a, (3) are the (dimensionless) expectation values and 
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pair correlations of two single spins {a, (3) to which it is assumed there is no experimental access: 



J ' Pcfi 

They are obtained from the actual 4x4 density matrix of the two spins a and /3: 

3 3 3 



(3.49) 
(3.50) 



Pa/S 



l(")0l('3) + ^5fc(a)ai")®l(« + ^l(")05/(/3)'^f^+E'^fc'("'/3)4"^^^r (3-51) 

fc=l 1=1 k,l=l 

with 1^"^ (1^'^)) the 2x2 identity matrix in the Hilbert space of spin a (/3). 

Out of the experimentally accessible quantities (3.47) and (3.481 we will construct a 4x4 density 
matrix of two virtual qubits which describes the collective properties of the two spin ensembles. The 
a priori justification for this method is: 

1. A general treatment of the problem between two large samples of spins is intractable because 
of the high dimensionality. 

2. We have a fully developed theory of entanglement for two-qubit systems. Therefore, this 
approach is a natural way to say something about the entanglement between two spin systems 
if only collective observables are measured. 

We first introduce the normalized (dimensionless) average subsystem expectation values (mag- 
netization per particle) and correlations: 



< 4 s »■<''> =Ji*- 



(3.52) 
(3.53) 



where s^,<^j' e [-1, 1]. These are the coefficients of the virtual density matrix: 




(3.54) 



with a denoting the first and b the second virtual collective qubit, associated with subsystems A 

7^ and (jf 



and B, respectively. Here, 1", l**, ct^ and trf are 2x2 identity and Pauli matrices for the collective 



qubits a and b. 

The question is whether the density matrix (3.54) is positive semi-definite, i.e. whether it is a 
physical state of two qubits. The answer is affirmative and the proof follows from the consideration 
of an equal-weight statistical mixture of one (virtual) qubit pair which can be in any of the 
states Pap- The density matrix of this mixture is the mixture of density matrices of all possible 
pairs (a, (3): 

Pmix = ^ E « • (3-55) 
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It can easily be seen that 

Pmix = Pab (3.56) 
for both are uniquely determined by the same expectations and correlations: 



{^f')p^^^{i^'i)p..=^^^St, (3.57) 
^r'-f),^.. - ko-o'^K. - - ^ ■ (3-58) 



Thus, Pah is indeed a density matrix. Note that without the normalizations as given in (3.52) and 



(3.531, the method would not work. 

Encapsulated in the following two propositions, we relate the entanglement properties of the 
virtual qubits to those of the spin samples. 

Proposition 1. For any entanglement measure E that is convex on the set of density matrices 
the entanglement of the virtual density matrix Eab = E(pab) is a lower bound for the average 
entanglement between all pairs Eap = Tla p ^{Pap)- 

Proof: This is an immediate consequence of the convexity of E: 



Remarks: First, the result holds for entanglement measures that are convex. In certain cases 
this is directly implied by the definition of the entanglement measure for mixed states, which 
involves a convex roof E{p) = min^.^y,. ^^pi i<^(|'0i)('0i|), where the minimization is taken over 
those probabilities pi and pure states |'0i) that realize the density matrix p = J2iPi\''Pi){''P'i\ ^^'^ 
E!{\ipi){tpi\) is the entanglement measure of the pure state Second, the proposition implies 

that if Eab > then at least for one pair (a,/3) we must have E{pap) > 0: 

/5ob entangled ^ 3(q;, /?): pa/j entangled . (3.60) 

Thus, a non-zero value of Eab is a sufficient condition for entanglement between the two samples A 
and B. Third, the maximal pairwise concurrence |42j for symmetric states is found to be 2/n and 
is achieved for the W-state [H]. It is conjectured that this remains valid also when the symmetry 
constraint is removed. This suggests Eab < Eaf} < 2/n, if concurrence is used as an entanglement 
measure. The existence of this upper bound can be seen as a consequence of the monogamy of 
entanglement. 

We refer to Eab £^s pairwise collective entanglement as it is determined solely by the expectation 
and correlation values of the collective spin observables. The question arises: Under what conditions 
is Eab equal to the average entanglement Eapl Identifying systems for which the equality holds, 
would allow feasible experimental determination of the entanglement distribution in large samples 
by observation of their macroscopic properties only. It can easily be seen that, if the state is 
symmetric under exchange of particles within each of the samples, one has Eab — Eap — E{pap) 
for every pair of particles (a, (3). In what follows, we identify an important class of systems for which 
Eab — Ea[3 though thc corresponding states need not be symmetric under exchange of particles. 
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Proposition 2. Consider a system (i) with gz{a) — gf for all a G A and gz{P) — gf for all 
(3 € B (translational invariance within the subsystems), (ii) with hxx{oi,P) — shyy{a, (3) with 
£ = const = +1 or —1 (all pairs are in absolute value equally correlated in the x and y-direction, 
where this correlation may be different in size for different pairs), (iii) with constant sign of the 
2;-correlations, i.e. sgn[hzz{ot, (3)] = — e for all (a,/3), and (iv) where all the remaining expectation 
values and correlations {gx, gy, hij with i ^ j) are zero. Non-vanishing average entanglement Eafj 
as measured by the negativity is equal to the pairwise collective entanglement Eab, if and only if 
the correlation functions hxx{oi,P) and hyy{a, (3) are each constant for all pairs and all pairs have 
a non-positive eigenvalue of their partial transposed density matrix. 

Proof: The negativity jl05l of a density matrix p is defined as 



(3.6f) 



where TrlpP"""! stands for the trace norm of the partially transposed density matrix pP""". Hence the 
negativity is equal to the modulus of the sum of the negative eigenvalues of pP""". 

It is important to stress that proposition 2 holds also for states that do not need to be totally 
symmetric, i.e. the hzz(oi,f3) may be different for different pairs of particles. In general, under the 
above symmetry, the state of the virtual qubit pair is of the form 



Pab 



( W++ 






(3.62) 



w_ \ 

77-1- 

H-i- u \ 

u—^j 

\ (lie). The state of an arbitrary pair (a, (3) of 



where u±±± = \ (l ± ± sf ± tfj, v± 
particles has a similar structure. For example, the state of an arbitrary pair of particles extracted 
from a spin chain with xxz Heisenberg interaction has such a form. 

Depending on the sign e of the zz-correlations, only one eigenvalue of and p^J , respectively, 
can be negative: 



Ma0 



1 

4 
1 



1 - ^{g^ + eg^r+AhUa,P) + e hzz{a, /?) 
l^^{s- + eslY+A{t%iY + et\ 
The corresponding negativities are given by 



lab 

^zz 



(3.63) 
(3.64) 



E{pap) = |min(0,pa/3)| 
Eab = |min(0, Vab) \ ■ 



(3.65) 
(3.66) 



One can express Vab as given by 



where 



Vab^ tJ' + 



(3.67) 
(3.68) 
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and 

^ ^ E„,^ V + ^ ff?)' + 4 hUa, /3) _ i + e + 4 {tf^Y . (3.69) 

The quantity A is the difference between the entanglement measures Eaf3 and Eab, i.e. Eab — 
Eap — A, for the case that Vab < and Eap = ;^ X^a /3 I™ii(Oj A'a/j)! = |min(0, /i)|. This is true, 
if and only if /^^^ < for all {a, (3), i.e. all pairs are either entangled or have eigenvalue zero. 
According to proposition 1, A is non-negative, i.e. 

miY<\Y. ,Vc^+^hUa,P). (3.70) 



Here we abbreviated c = gf + e — + £s^, where the latter equal sign is due to (3.521 



Inequality (3.701 becomes an equality, i.e. A = 0, if and only if hxx{oL^(3) is the same for all pairs 
(a,/3) such that — h^x- Therefore, the pairwise collective entanglement E^b equals the average 
entanglement i?Q/3, if and only if for all individual pairs jiap < and /ia;a;(Q;, /3) =ehyy{a,(3) = const 
for all pairs. □ 

We illustrate the method with some explicit examples: 

1. Dicke states 

We consider the Dicke state (generalized W-state) 

\N;k) ^ [^,y'^' Ps \0^1^) (3.71) 

N-k k 

with N > 2 spins, k excitations |1) and N — k non-excited spins |0), where < fc < A^. Ps is 



the symmetrization operator. Here, we identify the collective spin operators Si, eq. (3.45), which 
act on 2" dimensional Hilbert spaces, with spin operators Si acting only on the 2s + 1 = n + 1 
dimensions of the (symmetric) Dicke states. 

Within the system we consider two subsystems A and B each of size n. Because of the total 
symmetry of the state one has Eab = Eafi = E{pa[j) for any size n. Only for the cases where just a 
single spin (k = 0) or all spins are excited (fc = N), there is no entanglement between two arbitrary 
pairs or arbitrary sized blocks, respectively |92j . 

The reduced two-qubit density matrix has the form 

p^p{N- k) = d |00)(00| + e |11)(11| + 2/ |^+)(^+| , (3.72) 
where = TjdOl) + |10)) and 

iN-k)iN-k-l) ^ fc(fc-l) fc(jV-fc) 

N{N-1) ' A(A-l)' ^ N{N-1)' ^ ' 

with d + e + 2f = 1. The reduced density matrix pap is the same for all pairs of spins (a, (3), 
independent of their position and distance from each other. We need to calculate the expectation 
values of the spin operator and correlations between the two sites a and (3. They are independent 
of a and /3: 

9z = e-d, hxx^hyy^2f, h^^ = d + e-2f . (3.74) 
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All the others {gx,gy,hij with i ^ j) are zero. This is a Heisenberg xxz type situation, for which 
all our proofs above hold. The corresponding normalized collective values of two n-particle blocks 



A and B, (3.521 and (3.53), are 



= - ^^'^ = - y: 5. = 5. , (3.75) 



n n 



t^= = -2T^''--2ll 11^^==^^^^ (3-76) 

a^A l3eB 

i.e. the same as the actual values for two arbitrary individual spins. This means that for the W- 
state the averages are the values themselves, since there is no dependence at all on the position of 
the spins and therefore on the distance between them. Therefore, the collective block entanglement 
(for all block sizes n < ^) equals the entanglement between two (arbitrary) spins: 

Eab = E{p^p) . (3.77) 

The corresponding (possibly negative) eigenvalues of the partial transposed matrices are 

Vah = IJi<.p = \{.d+e)- \^{e-dY + Ap. (3.78) 

The entanglement 

Eab = Wab\ (3.79) 

is non-zero for all excitations 1 < fc < iV — 1 . Only for the cases where not a single spin or where all 
spins are excited, Eab = and there is no entanglement between two arbitrary pairs or arbitrary 
sized blocks, respectively. The global maximum of the entanglement is reached for fc = ^, which 
is the case where half of the spins are excited {N should be even here) and its value is 



^ Eab{k^ f ) - \ ^ , (3.80) 

for all block sizes n. It vanishes in the limit of an infinite chain N ^ oo. Figure [3?8| shows the 
block entanglement Eab ^ £^ function of the number of excitations k. 

2. Generalized singlet states 

The two subsystems A and B, each forming a spin s = f , are in a generalized singlet state: 

|V'> = -^|^ E (3-81) 

where \m) — |2s; s + m) denotes the eigenstates of the spin operator's z-component. The collective 
two-qubit coefficients are t^^ — — and the s^"^ and t^^ with i ^ j are all zero. The eigenvalue 
(3.641 of becomes Vab — ^^id the collective entanglement (negativity) is 



Eab^:^^]-. (3.82) 
It is non-zero for all sizes n of the subsystems and vanishes only in the limit n oo. This agrees 



with the statements in References [531 [71] that the generalized singlet state (3.811 can violate a 
Bell inequality no matter how large the spin is. 
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Figure 3.8: Entanglement Eab between two arbitrary sized spin subensembles A and B, where the system 



is in a Dicke state (3.71 \ of A'' particles with k excitations. The entanglement is computed as the negativity 
of two virtual qubits a and b, representing the subensembles. It is non-zero except for the cases where not 
a single spin or where all spins are excited. The maximum entanglement E^i^^ = | j^zr[ is reached for 
k = N/2. 



3. Generalized singlet state with an admixture of non-symmetric correlations 



Consider the state 



1 

= pIV'X^I + (1 -p) (g) 2 (|oi>a,/3=a(io| + |io)„,^=„(oi|) 



(3.83) 



with p e [0, 1]. This is a mixture of the generalized singlet state (3.81 1 and n perfectly z-correlated 
pairs {a, (3 — a). This state is not symmetric under particle exchange in the zz-correlations. The 
expectation values s^'^ and correlations tf.^ with i ^ j remain zero. The correlations tx 



til 



are reduced by a factor p compared to those of the state (3.811. The correlations in z-direction, 
however, are modified and read = —p ^^^^ — ^. Therefore, there is a critical number of particles 
He = \ ] , with [ . ] the ceiling function, beyond which there is no collective pairwise entanglement. 
Only for n < ric we have positive 



1 + p - n (1 - p) 
4n 



(3.84) 



Note that (3.831 is in accordance with proposition 2 and thus Eab — Eap- Figure 3.9 shows Eab as 



a function of the spin length s — and the mixing parameter p. Eab is non-zero in regions where 
p > and decreases inversely proportionally to s. 



3.2.2 Multi-partite entanglement 

Our method can be generalized to define multi-partite entanglement of M collective spins belonging 
to M separated samples Ai, ...,Am, each containing a large number of spins n (Figure 3.101. The 
collective spin of subsystem Ap is 

sf;^^ H-'y- (3-85) 
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Eab 




0.25 



Figure 3.9: Entanglement Eab between two collective spins in the generalized singlet state with an ad- 



mixture of non-symmetric noise (3.831 as a function of spin length s and proportion p of the singlet state 



in the mixture. The entanglement is non-zero for sufficiently large p and decreases inversely proportionally 
to s. 



A. 



Figure 3.10: Schematic of a multi-partite scenario. There are M separated spin subensembles A\, Am, 
each containing a large number of individual spins. 

Here p = 1, M and S {0, 2:= 1. y = 2, z = 3} and CTq"^^ — 1 denotes the 2x2 identity matrix. We 
suppose again that the experimenter measures the collective spin components (where all ip but one 
are zero) in each sample and their (higher-order) correlations (where more than one is unequal 
to zero) Tt:itr ^(^^^••^^^): 

^v.-.iM^^pf E - E ^n....M(«i,-,«Af), (3.86) 

where 

/^,,......(ai, ...,«./) ^ (3.87) 

denotes the actual spin expectation values and higher-order correlations in the physical system. 
These are 4*^ — 1 numbers. The virtual correlations are denoted as 

nM 

,a\...aM _ rpA\...AM 

''il-.-iM — (j^f^^M -^il-.-iM ■ 1,0.00; 

Note that Sq'' = n 1 is n times the 2x2 identity matrix. A normalized p-particle correlation t1^"'°'^ 
with p < M indeed scales with n^P as p (M ~p) is the number of subscripts unequal (equal) to 
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zero. The virtual 2*^x2*^ collective Af-qubit density matrix is given by 

41=0 ijf=0 



(3.89) 



with CTj"'' denoting the 2x2 Pauli or identity matrix of the p-th virtual collective qubit, associated 
with the n-particle subsystem Ap. 

We show that the virtual matrix Pai...aM ^ density matrix. Analogously to the two-subsystems 
case, we consider an equal-weight statistical mixture of one m-tuple of qubits (ai, q;m) which 
can be in any of the n*^ (mixed) states Pai...aM- The density matrix of this mixture is the mixture 
of density matrices of all possible m-tuples: 



-y 

We have 



Pmix — ,r y Pqi...qj, 
71 — ^CKi,...,a^ 



Pmix — Pai...aM 

as both give the same expectations and correlations: 

i.("Af)\ _ i~.ai ;^a.M\ _ j.ai...aM 

Hence, analogously to proposition 1, any convex multi-partite entanglement measure (e.g., M- 
way tangle) which is applied to the collective matrix, E{pa-^^...aM)^ gives a lower hound for the average 
multi-partite entanglement -Eai...QM- For states which are symmetric within the samples equality 
holds: E{pa-^...aM) — Eai...aM- For totally symmetric states this is also equal to the multi-particle 
entanglement of M extracted particles from the system. Since being a legitimate entanglement 
measure, the collective multi-partite entanglement E{pai...aM) obeys the usual constraints for en- 
tanglement sharing such as the CofFman-Kundu-Wootters inequality [24] . 

Importantly, in the examples considered above our entanglement measure scales at most with 
1 /n and vanishes in the limit of infinitely large subsystem sizes n. This is a generic property that 
follows from the commutation relation for normalized spins in this limit. Taking = ^ Si = 

Ea 4"'' one obtains 



hm_ [sx,Sy] ^ lim^ i — = . (3.90) 



h 

n — ►CJO^ n — *oo 2/1 

This is sometimes interpreted as suggesting that averaged collective observables, like the magneti- 
zation per particle, represent "macroscopic" or classical-like, properties of samples. Note, however, 
that for any n there are pairs between the subsystems so that the number of pairs multiplied 
by the pairwise collective entanglement can scale with n, showing the existence of entanglement for 
arbitrarily large n. 
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Chapter 4 



Mathematical undecidability and 
quantum randomness 

Summary: 

The mathematics of the early twentieth century was concerned with the question whether a 
complete and consistent set of axioms for all of mathematics is conceivable [H]. In 1931 Godel 
showed that this is fundamentally impossible |39]. In every consistent axiomatic system that is 
capable of expressing elementary arithmetic there are propositions which can neither be proved nor 
disproved within the system, i.e. they are "undecidable" . Via the halting problem Turing brought 
Godel's mathematical proof into the world of physical machines [9l]. Chaitin went even one step 
further and argued that mathematical undecidability is not bound to self-referential statements but 
arises whenever a proposition to be proved and the axioms contain together more information than 
the set of axioms itself [H] . 

Here we propose a new link between mathematical undecidability and quantum physics. We 
demonstrate that the states of elementary quantum systems are capable of encoding mathematical 
axioms. Quantum mechanics imposes an upper limit on how much information can be encoded in 
a quantum state |13lll02j . thus limiting the information content of the set of axioms. We show that 
quantum measurements are capable of revealing whether a given proposition is decidable or not 
within this set. This allows for an experimental test of mathematical undecidability by realizing in 
the laboratory the actual quantum states and operations required. We demonstrate experimentally 
both the encoding of axioms, using polarization states of photons, and that the decidability of 
propositions can be checked by performing suitable quantum measurements. We theoretically find 
and experimentally confirm that whenever a mathematical proposition is undecidable within the 
system of axioms encoded in the state, the measurement associated with the proposition gives 
random outcomes. Our results support the view that quantum randomness is irreducible ^20) and 
a manifestation of mathematical undecidability. 

Despite its overwhelming success, quantum physics is still heavily debated for its interpretation 
has remained to be unclear. One of the reasons is that, in contrast to all other theories, it is 
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lacking clear and unambiguous foundational principles. The link with pure mathematics might 
provide such a principle and can be seen as a novel approach towards a reconstruction of quantum 
theory [Si] . 

This chapter mainly bases on and also uses parts of Reference [7D] : 

• T. Paterek, R. Prevedel, J. Kofler, P. Klimck, M. Aspelmeyer, A. Zeilinger, and C. Brukner 
Mathematical undecidability and quantum randomness 
Submitted (2008). 



4.1 Logical complementarity and mathematical undecidability 
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4.1 Logical complementarity and mathematical undecidabil- 
ity 

We begin our argumentation following the idea of the information-theoretical formulation of Godel's 
theorem |21j : Given a set of axioms that contains a certain amount of information, it is impossible to 
deduce the truth value of a proposition which, together with the axioms, contains more information 
than the set of axioms itself. To give an example, consider Boolean functions of a single binary 
argument: 

xe{0,i} ^ y = /(.t) e {0, 1} (4.1) 

There are four such functions, ?/fe(fc = 0,l,2,3), shown in Figure [4T[ We shall discuss the following 
(binary) propositions about their properties: 

(A) "The value of /(O) is '0', i.e. /(O) = 0." 

(B) "The value of /(I) is '0', i.e. /(I) = 0." 

These two propositions are independent. Knowing the truth value of one of them does not allow 
to infer the truth value of the other. Ascribing truth values to both propositions requires two 
bits of information. If one postulates only proposition (A) to be true, i.e. if we choose (A) as an 
"axiom", then it is impossible to prove proposition (B) from (A). Having only axiom (A), i.e. only 
one bit of information, there is not enough information to know also the truth value of (B). Hence, 
proposition (B) is mathematically undecidahle within the system containing the single axiom (A). 
Another example of an undecidable proposition within the same axiomatic system is: 

(C) "The function is constant, i.e. /(O) = /(I)." 

Again, this statement cannot be proved or disproved from the axiom (A) alone because (C) is 
independent of (A) as it involves /(I). 

We refer to independent propositions to which one cannot simultaneously ascribe definite truth 
values — given a limited amount of information resources — as logically complementary propositions. 
Knowing the truth value of one of them (i.e. having the proposition itself or its negation as an 
axiom) precludes any knowledge about the others. For example, given the limitation to one bit of 
information, the three propositions (A), (B) and (C) are logically complementary. 

When the information content of the axioms and the number of independent propositions in- 
crease, more possibilities arise. Already the case of two bits as the information content is instructive. 
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Figure 4.1: The four Boolean functions y = fix) of a binary argument, i.e. fix) = 0, 1 with a:: = 0, 1. The 
diflerent functions are labeled by yt with fc = 0, 1, 2, 3. 
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Consider two independent Boolean functions fi{x) and /2(a;) of a binary argument. The two bits 
can be used to define truth values of properties of the individual functions or they can define joint 
features of the functions. An example of the first type is the following two-bit proposition: 

(D) "The value of /i(0) is '0', i.e. /i(0) = 0." 
"The value of /2(1) is '0', i.e. /2(1) = 0." 

An example of the second type is: 

(E) "The functions have the same values for argument '0', i.e. /i(0) = /2(0)-" 
"The functions have the same values for argument '1', i.e. /i(l) = /2(1)-" 

Both (D) and (E) consist of two elementary (binary) propositions. Their truth values are of the 
form of "vectors" with two components being the truth values of their elementary propositions. 
The propositions (D) and (E) are logically complementary. Given (E) as a two-bit axiom, all the 
individual function values remain undefined and thus one can determine neither of the two truth 
values of (D). 

The new aspect of multi-bit axioms is the existence of "partially" undecidable propositions, 
containing more than one elementary proposition only some of which are undecidable. An example 
of such a partially undecidable proposition within the system consisting of the two-bit axiom (D) 
is: 

(F) "The value of /i(0) is '0', i.e. /i(0) = 0." 
"The value of /2(0) is '0', i.e. /2(0) = 0." 

The first elementary proposition is the same as in (D) and thus it is definitely true. The impossibility 
to decide the second elementary proposition leads to (partial) undecidability of the proposition (F). 
In a similar way, proposition (F) is partially undecidable within the axiomatic system of (E). 



4.2 Physical implementation and quantum randomness 

The discussion so far was purely mathematical. We have described axiomatic systems (of limited 
information content) using properties of Boolean functions. Now we show that the undecidability 
of mathematical propositions can be tested in quantum experiments. To this end we introduce a 
physical "black box" whose internal configuration encodes Boolean functions. The black box hence 
forms a bridge between mathematics and physics. Quantum systems enter it and the properties of 
the functions are written onto the quantum states of the systems. Finally, measurements performed 
on the systems extract information about the properties of the configuration of the black box and 
thus about the properties of the functions. 

We begin with the simplest case of a spin-i particle (qubit) entering the black box (see Figure 



4.2 



for the experimental implementation) and a single bit-to-bit function f{x) implemented by the black 
box. Inside the black box two subsequent operations alter the state of the input spin. The first 
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|z+> = Polarizer 0° 
|x+>= Polarizer 45° 
|y+>= Polarizer 0°+ 
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<y| : QWP(S)45° 



Figure 4.2: The experimental setup. We use the polarization of single photons as information carriers of 
binary properties encoded by the configuration in the "black box" . The photons are initialized in a definite 
polarization state, belonging to one of the three complementary bases z,x,y (Preparation). The Boolean 
functions are realized within the black box by inserting half-wave plates (HWPs) which implement the 
product of Pauli operators a^'"' (Tz''^', eq. (4.21. The measurement apparatus consists of a quarter-wave 
plate (QWP) or HWP, followed by a polarizing beam-splitter (PBS) and two fibre-coupled single-photon 
detector modules. In this way, measurements in all complementary bases (2, x, y) can be realized. 



operation encodes the value of /(I) via application of ai , i.e. the Pauli z-operator taken to the 
power of /(I). The second operation encodes /(O) with ai'"'^\ i.e. the Pauli x-operator taken to the 
power of /(O). The total action of the black box is 

U = ai^'^ai^'\ (4.2) 

Consider the input spin to be in one of the eigenstates of the Pauli operator i™" o-™ ct" (with i 
the imaginary unit). The three particular choices (m,ri) = (0, 1), (1,0), or (1, 1) correspond to the 
three spin operators along orthogonal directions ct^ , cTj; , or CTj, = i dx , respectively. The input 
density matrix reads 

p= i[l + A™„i'""a;"a,:], (4.3) 
with Amn = ±1 and 1 the identity operator. It evolves under the action of the black box to 

C/p [/^ = i [1 + A,„„ (-l)"/(o)+™/(i) i™" a^] . (4.4) 

Depending on the value of n/(0) + to /(I) (in this chapter all sums are taken modulo 2), the 
state after the black box is either the same or orthogonal to the initial one. If one now performs 
a measurement in the basis of the initial state, i.e. the eigenbasis of the operator i™" ct™ ct", the 
outcome reveals the value of n/(0) -I- to /(I) and hence the measurement can be considered as 
checking the truth value of the proposition 



(G) "n/(0) + TO/(l) = 0.' 
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Figure 4.3: The hnk between mathematical undecidability and quantum randomness. 



Each of the three quantum complementary measurements az, o'^, or ay reveals the truth value 
of one of the independent propositions (A), (B), or (C), respectively. 

Independent of the initial state, we now identify the quantum measurement (m, n) with the 
question about the truth value of the corresponding mathematical proposition (G). The states 
that give a definite answer in the quantum experiment encode (G) or its negation as an axiom. 
A measurement quantum physically complementary to the one identified with (G) gives random 
results for those states, and the corresponding logically complementary proposition is undecidable 
within the one-bit axiom. 

We will show in general that whenever the proposition identified with the measurement is 
decidable (in the axiomatic system encoded by the state after the black box), the measurement 
outcome is definite, and whenever it is undecidable, the measurement outcome is random. This 
links mathematical undecidability and quantum randomness and allows to experimentally find out 



whether a proposition is decidable or not. The essence of this chapter is summarized in Figure 4.3 



A natural explanation of irreducible (objective) randomness of individual quantum outcomes 
arises. After leaving the black box the spin's quantum state encodes exactly one bit of information 
about f{x), namely the truth value of the proposition (G). In mathematical language, the system 
encodes a one-bit axiom. One bit is the maximum amount of information which can be carried by a 
single spin-i particle [43 } 1102] . Thus, if this bit represents the truth value of, say, proposition (A), 
therefore defining (A) or its negation as an axiom, the other two logically complementary propo- 
sitions, (B) and (C), are undecidable. This is because there is no information left for specifying 
their truth values. However, the spin can nevertheless be measured in the bases corresponding to 
(B) or (C) and — as in any measurement — will inevitably give an outcome, e.g. a click in a detec- 
tor. The clicks must not contain any information whatsoever about the truth of the undecidable 
proposition. Therefore, the individual quantum outcome must be random, reconciling mathemati- 
cal undecidability with the fact that a quantum system always gives an "answer" when "asked" in 
an experiment. This provides an intuitive understanding of quantum randomness, a key quantum 
feature, using mathematical reasoning. 

To find out whether a proposition is undecidable, it is necessary to repeat an experiment suffi- 
ciently many times, such that — even in the presence of unavoidable experimental imperfections — the 
two possible different outcomes occur significantly often. To reveal decidability, only one outcome 
has to appear again and again, but sufficiently many repetitions of the experiment are needed, 
because in practice the other outcome sometimes occurs due to experimental errors. 

We illustrate the link between mathematical undecidability and quantum randomness in an 
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Measurement in z basis 



/(O) = ' 



yo yi yi 
IVIeasurement in x basis 



yo yi y2 ya 
IVIeasurennent In y basis 



/(0)+/(l) = ? 



yo yi y2 ys 

Figure 4.4: We confirm tliat a measurement in tiie z basis gives the value of /(O) and similarly, measure- 
ments in X and y bases give the value of /(I) and /(O) + /(I), respectively. Shown are the experimental 
results for the cases where we prepared and measured the qubit in the same basis. Here, green (blue) bars 
represent counts in detector "0" (detector "1"), see Figure 4.2 The top graph, e.g., shows the measurement 



in the z basis: /(O) is '0' (green) for the functions yo and y\, and '1' (blue) for 3/2 and 2/3. 



experiment. The qubit is realized by two orthogonal polarization states of a single photon generated 
in the process of spontaneous parametric down-conversion [5S] (SPDC). This system is formally 
equivalent to a spin-i particle. The horizontal/vertical, -1-45°/— 45°, right /left circular polarization 
of the photon corresponds to eigenstates |z±), |a;±), and |y±) of the spin-i particle, respectively. 
We start by initializing the qubit in a definite polarization state by inserting a linear polarizer in the 
beam path (see Figure 4.2 1. The qubit then propagates through the black box in which the Boolean 
functions are encoded with the help of half- wave plates (HWP). The desired unitary transformation 
<^z {<^x) on the polarization states is implemented by a HWP at an angle 0° (45°) with respect to 
the z-basis. Subsequently, measurements of Gz, cr^, and ay, which test the truth value of a specific 
proposition, are performed as projective measurements in the corresponding polarization basis. 
Specifically, we use a polarizing beam-splitter (PBS) performing dz measurements whose output 
modes are fiber-coupled to single-photon detector modules and use wave plates in front of the 
PBS to rotate the measurement basis (see Figure 4.2 1. The truth value of the proposition now 
corresponds to photon detection in one of the two output modes of the PBS. 
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Figure 4.5: Using the setup of Figure 4.2 we input the qubit in a well-defined Pauli operator eigenstate 
\z+}, or \y+) into the black box, shown from top to bottom. The black box encodes two classical bits, 

/(O) and /(I), and the choice of the initial state determines which single bit, /(O), /(I), or /(O) + /(I), 
is read out. For every input state we measure in all three complementary bases, z [asking for /(O)], x 
[/(I)], and y [/(O) + /(I)], shown from left to right. The three measurements are related to three logically 
complementary questions (A), (B), (C) of the main text as indicated by the labels. This particular plot is 
the experimentally obtained data for the black box realizing the function yi . Similar results were obtained 
for the other black box configurations yo, 1/2, and j/3. Again, green (blue) bars represent counts in detector 
"0" (detector "1") giving the answer to the corresponding question. Each input state, after leaving the black 
box, reveals the truth value of one and only one of the propositions, i.e. it encodes a one-bit axiom. Given 
this axiom, the remaining two logically complementary propositions are undecidable. This undecidability 
is revealed by complete randomness of the outcomes in the other two measurement bases. 



First, we confirm that complementary quantum measurements indeed reveal truth values of 
logically complementary propositions. To this aim we prepare the system in a state belonging to 
the basis in which we finally measure. The results are presented in Figure [4^ 

Next, for each of the three choices of the initial state, we ask all three logically complementary 



questions by measuring in all three different complementary bases. In Figure 4.5 we plot the count 
rate of photons measured after the PBS for a fixed configuration in the black box that encodes 



function yi. Similar results are obtained for other black box configurations (not shown). Figure 4.5 



shows that for every input state one and only one question has a definite answer, as already indicated 



in Figure 4.4 This is the axiom encoded in the system after it leaves the black box. The remaining 
propositions are undecidable given that axiom, and the corresponding measurement outcomes are 
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completely random, i.e. evenly distributed. 

4.3 Generalization to many qubits 

We now generalize the above reasoning using multiple qubits. Consider a black box with N input 
and N output ports, one for each qubit. It encodes N Boolean functions fj{x) numbered by 
j — 1 , . . . , iV by applying the operation 

Un = (jt^°^ H'^^^ H"^^'' ■ (4.5) 

The initial A'^-qubit state is chosen to be a particular one of the 2^ eigenstates of independent and 
mutually commuting tensor products of Pauli operators, numbered hy p = I, N: 

^ jmi(p)«i(p) ^mi(p) ^ . . . ^ jm«(p)n„(p) ^m„(p) ~«„(p) ^ 

with mj{p),nj{p) G {0,1}. A broad family of such states is the family of graph states [7S]- As 
before, the qubits propagate through the black box. After leaving it, their state encodes the truth 
values of the following N independent binary propositions (negating the false propositions, one has 
N true ones which serve as axioms): 

(Hp) "Ef=iK(p) m + m,{p) /,(!)] = 0." 

with p = 1, . . . ,N. Mathematically, N binary arguments being the N truth values of the (Hp) 
allow to construct 2 different Boolean functions. The statements about the values of these func- 
tions form all possible decidable (not independent) binary propositions. In suitable measurements 
quantum mechanics provides a way to test whether a given proposition is decidable or not. If one 
measures the operator 

e = i"^''^ CT^i fffi (g) • • • (g) i"""" a^" af" , (4.7) 
with aj,f3j G {0, 1}, one tests whether the proposition 

(J) "E,ti[/3j/.(0) + a,7,(l)]=0." 

is decidable or not. In case can be written as a product Cl^^ ■ ■ ■ Cl'^ with kp G {0, 1}, it commutes 
with all the fip's and consequently the measurement of Q has a definite outcome. On the other hand, 
proposition (J) is decidable within the set of axioms (Hp). Its truth value can be (logically) derived 
from the axioms in the sense that the sum in (J), X^^il/^i fji^) ^ linear combination 

of the sums in (Hp) with binary coefhcients kp. There are 2^ such linear combinations and therefore 
2^ such Q's. Since in general 4^ different Q's exist, the remaining 4^ — 2^ — 2^(2^ — 1) operators 
are linked with undecidable propositions and their measurement outcomes are random. Hence, there 
are much more undecidable propositions of the form (J) than decidable ones. The ratio between 
their numbers increases exponentially with the number of qubits, i.e. - — ^^jv~"'"'* = 0{2^). 

Interestingly, in case of (J) being decidable, its truth value imposed by (classical) logic is not 
necessarily the same as found in the quantum measurement. This can be demonstrated for three 
qubits initially in the Greenberger-Horne-Zeilinger (GHZ) state |40J 

|GHZ) = ^ i\z+), \z+), \z+), + \z-), |z-)3) , (4.8) 
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where e.g. \zzt)-^ denotes the eigenstate with the eigenvalue ±1 of az for the first qubit. We choose 
as axioms the propositions 

(KJ "/l(0) + /l(l) + /2(0)+/2(l) +/3(1) = 1." 

(K^) "/l(0) + /l(l) +/2(l) + /3(0) + /3(l) = l." 

(K3) " /l(l) + /2(0) + /2(l) + /3(0) + /3(l) = l." 

linked with the operators ay ® CTy ® cTx, cfy ® (t^ ® dy, and (ix ® (^y ® (^y, respectively. One can 
logically derive from (Ki) to (K3) the true proposition 

(L) "/l(l) + /2(l)+/3(l) = l." 

On the other hand, the proposition (L) is identified with the measurement of (Tx ® (Jx ® (Tx- But 
the result imposed by quantum mechanics corresponds to the negation of (L), namely: "/i(l) + 
/2(1) + /3(1) = 0." This is the heart of the experimentally confirmed GHZ argument against local 
realism (JUJ [BH EH] . In the (standard logical) derivation of (L) the individual function values are 
well defined and the same independently of the axiom in which they appear. Since this is equivalent 
to the assumptions of local realism, the truth values of decidable propositions found in quantum 
experiments do not necessarily have to be the same as the ones derived by logic. However, this does 
not change the connection between decidability (undecidability) of propositions and definiteness 
(randomness) of measurement outcomes. 

The question arises whether one can construct a classical device to reveal the undecidability 
of propositions. This is possible, provided one uses 2N classical bits to simulate quantum com- 
plementary measurements on N qubits [90]. Since such a device satisfies local realism, it gives 
the truth values of decidable propositions according to classical logic. On the level of elementary 
physical systems, however, the world is known to be quantum. It is intriguing that nature supplies 
us with quantum systems which can reveal decidability but cannot be used to learn the classical 
truth values. 

In order to illustrate the concept of partial undecidability as introduced above, we performed 
a two-qubit experiment. The two bits of proposition (E) described above correspond to the set 
of independent commuting operators f2i = az ® (Jz and ^2 = ® S'x- The common eigen- 
basis of these operators is spanned by the maximally entangled Bell states (basis 6e): \^^} = 
^ {\z+)^ |z+)2 ± \z-)^ \z-)2), 1^''^) = ^ {\z+)i \z~)2 ± \z-)i \z+)2)- Thus, after the black box 
the four Bell states encode the four possible truth values of the elementary propositions in (E) and 
a so-called Bell State Analyzer (i.e. an apparatus that measures in the Bell basis) reveals these 
values. In the same way, the truth values of the elementary propositions in (F) are encoded in the 
eigenstates of local az bases, i.e. by the four states \zzt)^ l-z±)2 (basis bp). Finally, the elementary 
propositions in (D) are linked with the product states Izi)^^ |3;±)2 (basis 6d). 

The setup that was employed to generate maximally entangled Bell states via type-II SPDC is 



shown in Figure 4.6 We also prepared separable two-qubit states by inserting polarizers into the 
optical paths just before the black box. The encoding of functions within the black box as well 
as the polarization measurements on the individual qubits were done similarly to the single-qubit 
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(a) (b) 

"Cond 1" 




"Cond 2" 

Figure 4.6: The setup in (a) allows for separable two-qubit measurements on product as well as on en- 
tangled input states. Joint two-qubit measurements require a so-called Bell State Analyzer (BSA), whose 
experimental realization is depicted in (b). In both cases, an ultra-violet laser pulse passes through a 
non-linear crystal (BBO) to produce polarization-entangled photon pairs in the process of spontaneous 
parametric down-conversion. Compensators (Comp), made up of half- wave plates (HWP) and BBO crys- 
tals, are used to counter walk-off effects in the down-conversion crystal. They are set such that \^^) states 
are emitted. In (a), local measurements in all complementary bases are performed with the help of HWP, 
QWP and PBSs. Measurements in the Bell basis, as depicted in (b), require an ancillary entangled Bell 
state I'E'^). Therefore the UV-laser passes twice through the crystal, emitting the input Bell state in the 
forward direction (modes a & b), and the ancilla pair in the backward direction (modes c & d). Coherently 
combining these photons on PBSs, allows the identification of all four Bell states whenever there is one 
photon in each output mode of the PBSs. Then, conditioned on the detection of one photon with 4-45° 
polarization in both detectors, "Cond 1" and "Cond 2" , the Bell state can be identified by analyzing the 
remaining two photons in the -f45°/— 45° basis (see Reference |95] for details). 



case. However, this only allows for independent single-qubit measurements. Joint measurements 
on both qubits in the Bell basis require conditional operations on two qubits. The heart of such 



a Bell State Analyzer [SH EH [^5] (BSA) is a controUed-NOT gate [i [TSJ [Ml [31] • In Figure 4.6 



we depict the experimental setup implementing the BSA. Details of its working can be found in 
Reference [55] . 

We start our investigation by confirming that the truth values of the (elementary) propositions in 
(E), (F) and (D) are revealed by measurements performed in the bases 6e, bp and bu, respectively. 
As in the single-qubit case, preparation and measurement are in the same basis. The obtained 
results (not shown) arc similar to Figure 4.4 Moreover, one can initially prepare a Bell state and 
measure it in the bases 6e, bp and 6d- As can be seen in the left plot of Figure 4.7 measurements 
in the Bell basis, 6e, prove that the entangled state indeed encodes joint properties of the functions 
/i(x) and ,f2{x), i.e. information about (E). Measurements in other bases can then be interpreted in 
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Figure 4.7: In this two-qubit experiment a l"!"^) Bell state is measured in three different bases (the 
Bell basis as well as l^i)2 ^'^'i l^i)2' shown from left to right) using the setups depicted in 

Figure [4^ Plotted are the count rates associated with different detector combinations for the black box 
encoding function y2 on both photons. Label "00" corresponds to a coincidence event in detectors "Oa" 
and "Ob" , and similarly for the other labels. The first (second) label gives the answer to the upper (lower) 
question. The measurements performed in the Bell basis show that the two qubits encode proposition (E) 
of the main text. The data in the middle plot reveal partial undecidability of proposition (F), given (E) as 
an axiom, as indicated by the random outcomes in two out of four detector combinations. In contrast, the 
right plot presents the data corresponding to the fully undecidable proposition (D), where the outcomes 
are completely random. Similar results for propositions (E), (F) and (D) were obtained for other black box 
encodings. The reason for the fact that the Bell state in the left plot is not identified with unit fidelity 
stems from imperfections in the experimental setup. Unequal detector efficiencies explain the small bias in 
the right plot. 



terms of partial and full undecidability. Proposition (D) is fully undecidable given (E) as an axiom 
which is encoded by !<&''') after the black box. This can be seen from the right part in Figure 



4.7 



which all four measurement outcomes occur with equal probability. On the other hand, proposition 
(F) is partially undecidable. This is experimentally revealed by the count distribution of the middle 



part in Figure 4.7 The partial undecidability is uncovered by the randomness of the two occurring 
outcomes, while the other two outcomes do not appear. 

In conclusion, we have demonstrated how mathematical axioms can be encoded in quantum 
states and how decidability of mathematical propositions can be verified in quantum measure- 
ments. This extends the concept of "experimental mathematics" [22] to a new domain and sheds 
new light on the (mathematical) origin of quantum randomness. We have performed an experiment 
which showed for the first time the link between mathematical undecidability and quantum ran- 
domness. There, we found that whenever the quantum system is measured in a basis associated to 
an undecidable proposition, it gives random outcomes. Our results can be extended in many ways. 

An interesting avenue for further research arises from the observation that in mathematical 
logic one usually makes a difference only between statements that arc cither decidable or undecid- 
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able within a formal system of axioms. Quantum complementarity implies that losing certainty 
about one of the propositions is followed by a corresponding gain of certainty in other, logically 
complementary, propositions |16l 117] . This opens up the possibility to quantify the amount of 
undecidability of a proposition from "impossible" to "necessary" , and to describe the continuous 
trade- off in gaining and losing knowledge about them. 
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This final part of the present work offers the room to point out some possibihties for further research, 
including not only conservative but also speculative thoughts. 

Macroscopic realism and the quantum-to-classical transition: 

• We have put forward an approach to the quantum-to-classical transition, resting solely on 
the idea that coarse-grained measurements give rise to the emergence of macroscopic realism. 
However, to offer a complementary approach to the decoherence program that also applies to 
isolated systems, it is necessary to generalize the formalism to systems other than spins. E.g., 
for a particle in conventional phase space with position and momentum and an arbitrary 
time evolution there is no quantum number that can be made larger and larger (like spin 
length). But Reference |74j shows how a small blurring of a particle's Wigner function makes 
it positive (by convoluting it with the Wigner function of a coherent state). Again, it should 
be possible to model coarse-grained measurements in such a way that they are non-invasive at 
the level of a positive probability distribution. The latter represents an ensemble of classical 
objects with definite positions and momenta and must be sufficient to compute probabilities 
for coarse-grained outcomes. It can be seen already now that the interesting issue of non- 
classical Hamiltonians will arise again. 

• We have found that coarse-grained spin measurements are best modeled with a positive oper- 
ator value measure because coarse-grained von Neumann measurements allow to distinguish 
microstates at two sides of a slot border. Still, under all circumstances it is unavoidable that 
a quantum measurement always introduces a slight disturbance of the state even on the level 
of the Q-distribution. Under ideal experimental conditions and sufficiently many runs, we 
should be able to see tiny deviations from classical physics even for classical Hamiltonians, 
macroscopic systems, and coarse-grained measurements. 

• We have demonstrated that every non-trivial Hamiltonian allows to violate the Leggett-Garg 
inequality as long as sharp quantum measurements can be performed. This might be seen 
as related to the quantum Zeno effect [66]: Quantum measurements not only are capable of 
freezing the time evolution of a system but they also may influence it in such a way that it 
cannot be described classically anymore. Note that a perfect quantum Zeno effect itself can 
easily be simulated classically just by using no time evolution at all. 

• We have seen that non-classical Hamiltonians can lead to a violation of macrorealism even 
under the restriction of coarse-grained measurements. Moreover, even environmental deco- 
herence, though leading to macrorealism, does not resolve the problem that in principle no 
classical laws of motion are able to describe such time evolutions. We have tried to argue 



96 



Conclusions and outlook 



why non-classical Hamiltonians do not appear in nature. One avenue of research would be 
to investigate the relation between the interactions in our measurement apparatuses and the 
ones governing the time evolutions of systems. The Hilbert space has no notion of closeness 
or distance of orthogonal states. We expect that under coarse-grained measurements which 
bunch together states by whatever notion of closeness or distance quantumness can be seen 
only if the system Hamiltonian and the apparatus Hamiltonian connect states differently. 

• We may even conjecture that the three- dimensionality of our real space is a necessity under 
the postulate that an intersubjective classical world arises out of the quantum realm merely 
due to coarse-graining. This is supported by the fact that the Lie algebras SU(2) and S0(3) 
are isomorphic. 

Quantum randomness and mathematical undecidability: 

• In the spirit of the approach in Reference |17j . one important question is whether it is possible 
to derive Malus's law in the framework of propositions, i.e. a very particular kind of contin- 
uous trade-off between gaining and losing certainty about truth values like in the quantum 
case. It might be that to this end the necessary number of independent binary functions 
immediately becomes very large and is related to the number of spatial directions that can 
be experimentally distinguished, or that one cannot use binary functions any longer. 

• Our link between quantum physics and pure mathematics offers an intuitive understanding 
of randomness in which the question of decidability is more fundamental than about truth. 
Which other ingredients beyond mathematical undecidability are necessary to fully reconstruct 
the quantum formalism of state vectors in a Hilbert space? Is it possible to derive from this link 
any other physical phenomena that cannot be derived or understood from the conventional 
quantum postulates? 

• Godel's incompleteness theorem is known to be equivalent to Turing's halting problem. The 
latter can be seen as a consequence of the fact that there are more real numbers than natural 
numbers, or equivalently, not uncountable many Turing machines for countably many possible 
inputs. Chaitin's version rests on the intuitive observation that a theorem cannot be derived 
from axioms if the information content of theorem plus axioms is larger than the information 
content of the axioms itself. Thus, in both cases it is a mismatch of "resources for answers" 
and "amount of possible questions" leading to undecidability — just like in quantum mechan- 
ics where randomness is a consequence of the fact that we can measure more things than 
the state can definitely define. Formalizing this link between Chaitin's and Godel's original 
version would strengthen the connection between quantum randomness and mathematical 
undecidability. 
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